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ABSTRACT 

This chapter reviews theoretical work on the stellar dynamics of galaxy disks. All the known collective 
global instabilities are identified, and their mechanisms described in terms of local wave mechanics. A 
detailed discussion of warps and other bending waves is also given. The structure of bars in galaxies, 
and their effect on galaxy evolution, is now reasonably well understood, but there is still no convincing 
explanation for their origin and frequency. Spiral patterns have long presented a special challenge, and 
ideas and recent developments are reviewed. Other topics include scattering of disk stars and the survival 
of thin disks. 

Subject headings: galaxies: evolution - galaxies: halos - galaxies: kinematics and dynamics - galaxies: 
spiral 
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1. INTRODUCTION 

A significant fraction of the stars in the universe resides 
in the rotationally supported disks of galaxies. Disks are 
mostly thin and flat, but the disk is often warped away 
from its principal plane in the outer parts. Disk galax- 
ies usually manifest spiral patterns, and rather more than 
half host bars. Most, but not all, disk galaxies have a cen- 
tral bulge, perhaps also a thick stellar disk, and generally 
a small fraction of the stars resides in a quasi-spherical 
stellar halo, while the central attraction at large distances 
from the center is dominated by a dark halo. The mate- 
rial in most disks overwhelmingly orbits in a single sense, 
although a small fraction of galaxies have been found to 
host substantial counter-rotating components. 

This chapter is primarily concerned with the dynam- 
ics of rotationally supported disks of stars. Stellar disks 
are believed to have formed over time from gas that had 
previously settled into centrifugal balance in the gravita- 
tional well of the galaxy, and the process of star formation 
continues to the present day in most disk galaxies. While 
stars are the dominant dynamical component today, the 
small gas fraction (usually < 10% by mass) can still play 
an important dynamical role in some contexts. 

Disk dynamics is a rich topic for two principal reasons: 
(a) the organized orbital motion facilitates gravitationally- 
driven collective behavior and (b) outward transfer of an- 
gular momentum extracts energy from the potential well. 
Space limitations preclude a detailed development and this 
review will mostly be confined to a summary of the prin- 
cipal results and open issues. The derivations of the prin- 
cipal formulae can be fou nd in the excellent textbook by 
iBinnev fc Tremaind (|2008l . hereafter BT08). Furthermore, 
no attempt is made to cite every paper that relates to a 
topic. 

The distribution of gas within the Milky Way, the prop- 
erties of the Galactic bulge, and other closely related topics 
are described elsewhere in this volume. The distributions 
of light and mass within galaxies and our current under- 
standing of the processes that lead to the formation of 
galaxies are described in volume 6. 



2. PRELIMINARIES 

2.1. Relaxation Rate 

Because stellar disks contain many stars, the attraction 
from individual nearby stars is negligible in comparison 
with the aggregated gravitational field of distant matter. 
The Appendix explains how the usual rough calculation 
to support this assertion, which was derived with quasi- 
spherical systems supported by random motion in mind 
(e.g. BT08 pp. 35-38), must be revised for disks. Three 
factors all conspire to reduce the relaxation time in disks 
by several orders of magnitude below the traditional esti- 
mate fea. lA4p . although it remains many dynamical times. 

Fig. [H shows the non-smooth distribution of stellar ve- 
locities of > 14 000 F & G dwarf stars in the vicinity of the 
Sun, as fou nd in the Geneva-Copenhagen Survev ( Nord- 
strom et al. l2004l : iHolmberg et a;.ll2009l . hereafter GCS). 
The determination of the radial velocities has confirmed 
the substructure that was first identified by Dehncn (1998) 
from a clever analysis of the HIPPARCOS data without the 
radial velocities (The implications of this Figure are dis- 
cussed more fully in jjlOl ) As collisional relaxation would 
erase substructure, this distribution provides a direct illus- 
tration of the coUisionless nature of the solar neighborhood 
(unless the s ubstructure is being recreated rapidlv. e. a. De 
Simone et al. l2004f ). 

The usual first approximation that stars move in a smooth 
gravitational potential well therefore seems adequate. 

2.2. Mathematical Formulation 

This approximation immediately removes the need to 
distinguish stars by their masses. A stellar system can 
therefore be described by a distribution function (DF), 
f{x, v,t) that specifies the stellar density in a 6D phase 
space of position x and velocity v at a particular time 
t. Since masses are unimportant, it is simplest conceptu- 
ally to think of the stars being broken into infinitesimal 
fragments so that discreteness is never an issue. 

With this definition, the mass density at any point is 
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Fig. 1. — The velocity distribution, in Galactic coordinates, of 
stars near the Sun, as given in the GCS. U is the velocity of the 
star toward the Galactic center and V is the component in the di- 
rection of Galactic rotation, both measured relative to the Sun. The 
intersection of the vertical and horizontal lin es shows the local stan- 
dard of rest estimated bv lAumer fc BinnevI feOOQ.') . 



which in turn is related to the gravitational potential, <J>, 
through Poisson's equation 

V^^ix,t) ^ A7rGp{x,t). (2) 

This is, of course, just the potential from the stellar com- 
ponent described in /; the total potential includes con- 
tributions from dark matter, gas, external perturbations, 
etc. Finally, the evolution of the DF is governed by the 
collisionless Boltzmann equation (BT08 eq. 4.11): 



df ^ df.df 

h V ■ h V ■ 

dt dx dv 



0, 



(3) 



where the acceleration is the negative gradient of the 
smooth total potential: i» = — V<i>tot- The time evolution 
of a stellar system is completely described by the solution 
to these three coupled equations. Note that collisionless 
systems have no equation of state that relates the density 
to quantities such as pressure. 

The most successful way to obtain global solutions to 
these coupled equations is through A^-body simulation. 
The particles in a simulation are a representative sam- 
ple of points in phase space whose motion is advanced in 
time in the gravitational field. At each step, the field is 
determined from a smoothed estimate of the density dis- 
tribution derived from the instantaneous positions of the 
particles themselves This rough and ready approach is 
powerful, but simulations have limitations caused by noise 
from the finite number of particles, and bias caused by the 
smoothed density and approximate solution for the field, 
and other possible artifacts. 

Understanding the results from simulations, or even 
knowing when they can be trusted, requires dynamical in- 
sight that can be obtained only from analytic treatments. 
This chapter therefore stresses how the basic theory of 
stellar disks inter-relates with well designed, idealized sim- 
ulations to advance our understanding of these complex 
systems. 

^ i.e. a simulation solves eqs. (l)-(3) by the method of characteristics. 



Fig. 2. — An orbit of a star is generally a non-closing rosette 
(upper panel). The lower panel shows one radial oscillation, from 
pericenter to pericenter say, during which time the star moves in 
azimuth through the angle A</>. 



2.3. Orbits 

When orbital deflections by mass clumps can be ne- 
glected ( ji2.ip . the stars in a disk move in a smooth grav- 
itational potential. The effects of mass clumps and other 
non-uniformities are considered in §10. It is simplest to 
first discuss motion of stars in the disk mid-plane, before 
considering full 3D motion. 

The orbit of a star of specific energy E and specific 
angular momentum in the mid-plane of an axisym- 
metric potential is, in general, a non-closing rosette, as 
shown in Fig. [2] The motion can be viewed as a retro- 
grade epicycle about a guiding center that itself moves 
at a constant rate around a circle of radius which is 
the radius of a circular orbit having the same L^- In one 
complete radial period, the star advances in azimuth 
through an angle A^, as drawn in the lower panel. In fact, 
TT < A(j) < 27r in most reasonable gravitational potentials; 
specifically Acf) = y/2TT for small epicycles in a flat rotation 
curve. 

These periods can be used to define two angular fre- 
quencies for the orbit: fl^ = lS.(j)/TR, which is the angular 
rate of motion of the guiding center, and ^Ir — 2tt /tj^. In 
the limit of the radial oscillation amplitude, a — 0, these 
frequencies tend to VL^ — > f2, the angular frequency of a 
circular orbit at R = Rg, and flu k, the epicylcic 
frequency, defined through 
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2.4. Resonances 
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Fig. 3. — The Lindblad diagram for a disk galaxy model. Circular 
orbits lie along the full-drawn curve and eccentric orbits fill the 
region above it. Angular momentum and energy exchanges between 
a steadily rotating disturbance and particles move them along lines 
of slope Qp as shown. The dotted and dashed lines are the loci of 
resonances for an m = 2 perturbation of arbitrary pattern speed. 

If the potential includes an infinitesimal non-axisym- 
metric perturbation having m-fold rotational symmetry 
and which turns at the angular rate ilp, the pattern 
speed, stars in the disk encounter wave crests at the 
Doppler-shifted frequency m|r2p — il^j. Resonances arise 
when 

m{flp - fl^) = IQr, (5) 

where simple orbit perturbation theory breaks down for 
steady potential perturbations. At corotation (CR), 
where I = 0, the guiding center of the star's orbit has the 
same angular frequency as the wave. At the inner Lind- 
blad resonance (ILR), where I — —1, or at the outer 
Lindblad resonance (OLR), where I = +1, the guiding 
center respectively overtakes, or is overtaken by, the wave 
at the star's unforced radial frequency. Other resonances 
arise for larger but these three are the most important. 
Note that resonant orbits close after m radial oscillations 
and / turns about the galactic center in a frame that ro- 
tates at angular rate flp. 

The solid curve in Fig. [31 the so-called Lindblad dia- 
gram, marks the locus of circular orbits for stars in an 
axisymmetric disk having a flat rotation curve. Stars hav- 
ing more energy for their angular momentum pursue non- 
circular orbits. The broken curves show the loci of the 
principal resonances for some bi-symmetric wave having 
an arbitrarily chosen pattern speed; thus resonant orbits, 
that close in the rotating frame, persist even to high ec- 
centricities. 

2.5. Wave-Particle Scattering 

The arrows in Figure [3] indicate how these classical inte- 
grals are changed for stars that are scattered by a steadily 
rotating mild potential perturbation. Since Jacobi's in- 
variant, 

Ij = E~ rtpL, (6) 

(BT08 eq. 3.112), is conserved in axes that rotate with the 
perturbation, the slope of all scattering vectors in Fig. |3] 



is AE/AL^ = n„. 

Lynden-Beil & Kalnaj^ ()1972[ ) showed that a steadily 
rotating potential perturbation causes lasting changes to 
the integrals only at resonances. Since flc = at CR, 
scattering vectors are parallel to the circular orbit curve 
at this resonance, where angular momentum changes do 
not alter the energy of random motion. Outward transfer 
of angular momentum involving exchanges at the Lindblad 
resonances, on the other hand, does move stars away from 
the circular orbit curve, and enables energy to be extracted 
from the potential well and converted to increased random 
motion. 

Notice also from Fig. [3] that the direction of the scatter- 
ing vector closely follows the resonant locus (dotted curve) 
at the ILR only. Thus, when stars are scattered at this 
resonance, they stay on resonance as they gain random 
energy, allowing strong scatterings to occur. The oppo- 
site case arises at the OLR, where a small gain of angular 
momentum moves a star off resonance. 

2.6. Vertical Motion 

Stars also oscillate vertically about the disk mid-plane. 
If the vertical excursions are not large, and the in-plane 
motion is nearly circular, the vertical motion is harmonic 
with the characteristic frequency v = (9^$/9z^)^/^, which 
principally depends on the density of matter in the disk 
mid-plane at Rg. As v is significantly larger than k, the 
frequency of radial oscillation, the two motions can be con- 
sidered to be decoupled for small amplitude excursions in 
both coordinates. In practice, the vertical excursions of 
most stars are large enough to take them to distances 
from the mid-plane where the potential departs signifi- 
cantly from h armonic (i.e. for \z \ > 20 pc in the solar 
neighborhood, iHolmberg k Flvnnll2004( ) . 

The principal vertical resonances with a rotating den- 
sity wave occur where m(f2p — il^) = ±z/, which are at 
radii farther from corotation than are the in-plane Lind- 
blad resonances. Since spiral density waves are believed 
to have substantial amplitude only between the Lindblad 
resonances, the vertical resonances do not affect density 
waves, and anyway occur where the forcing amplitude will 
be small. 

The situation is more complicated when the star's ra- 
dial amplitude is large enough that v varies significantly 
between peri- and apo-galacticon. Also, the vertical peri- 
ods of stars lengthen for larger vertical excursions, allowing 
vertical resonances to become important (see ^A^. 

3. LOCAL THEORY OF DENSITY WAVES 

3.1. Plane Waves in a Thin Mass Sheet 

Since solutions of Poisson's equation can be obtained 
analytically only for the simplest mass distributions, the- 
oretical treatments necessarily make quite drastic simpli- 
fying assumptions. Local theory, summarized in this sec- 
tion, is built around a WKB-type potential solution for 
density variations in the form of short wavelength plane 
waves in a locally uniform, razor thin disk. In this case, 
a plane- wave disturbance Eie**'^ to the surface density in 
the (a;, y)-plane at some instant gives rise to the potential 

$i(x) = -^Eie^'^^e-I'^^l, (7) 
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(c/. eq. 6.31 of BT08). This relation strictly applies only to 
straight waves of infinite extent. However, the sinusoidal 
density variations cause the field of distant waves to cancel 
quickly, and the formula is reasonably accurate near the 
center of a short wave packet. Note that eq. ([7]) does not 
depend on the inclination of the wavefronts to the radial 
direction, but the curvature of the wavefronts must also 
be negligible, which requires |fci?| ^ 1, with R being the 
distance from the disk center. 

3.2. Dispersion Relations 

A dispersion relation gives the relationship between the 
wavenumber k and frequency uj of self-consistent waves. 
Assuming (a) the above relation between density and po- 
tential for a razor-thin disk and (b) the wave is of small am- 
plitude, so that second order terms are negligible, the local 
dispersion relation for short axisymmetric density waves 
may be written 

uj^ = -27rG^k\T{s,x) (8) 

()Kalnaislll965[ ). where S is the local undisturbed surface 
density and the complicated reduction factor J^, with 
two dimensionless arguments, is explained below. Lin & 
Shu (119661 ) independentlv derived a generalized relation 
for tightly-wrapped spiral waves, which required two addi- 
tional assumptions: (c) that the wave vector k is approx- 
imately radial, and (d) that the disturbance is not close 
to any of the principal resonances. Their better-known re- 
lation simply replaces the frequency of the wave with the 
Doppler-shifted forcing frequency experienced by the orbit 
guiding center: 

{tu -mQ)^ ^ k"^ ~27TGi:\k\Tis,x) (9) 

(see also BT08 §6.2.2), where now ui = mflp. 

The reduction factor, is always positive and is unity 
when the disk has no random motion. Since all the fac- 
tors in the self-gravity term on the RHS are intrinsically 
positive, eqs. ([5]) & (0) say that self-gravity enables a sup- 
porting response from the stars at a frequency that is lower 
than their free epicycle frequency k. 

When the disk has random motion, the vigor of the 
stellar response depends on two factors: the ratio of the 
forcing frequency experienced by a star to its natural fre- 
quency, s = \u} — mfl\/K, and the ratio of the typical sizes 
of the stellar epicycles (oc {v\)^^'^ / k) to the wavelength of 
the wave (a l^l"""^). Thus for a Gaussian distribution of ra- 
dial velocities, is a function of both s and x = I ■ 
Clearly, when x is large, the unforced epicyclic amplitude 
of most stars is larger than the radial wavelength, and 
the weakened supporting response arises mainly from the 
small fraction of stars near the center of the velocity dis- 
tribution. 

The dispersion relation for barotropic fluid (gaseous) 
diskfl is similar (BT08 eq. 6.55) 

(w -mf7)2 = - 27rGS|fc| -f i;^fc2^ (10) 

where Va is the sound speed in the gas. Both dispersion 
relations are shown graphically in Fig. 21 below. 

^ IGoldreich fc Lynden-Belll l ll965al) derived the vertically-integrated 
equations for "2D pressure" in a gas disk. 



3.3. Axisymmetric Stability 

Toomre (|1964f ) determined the condition for marginal 
stability of short axisymmetric waves by solving for = 
in eq. (|H]). Since all the factors in the self- gravity term are 
intrinsically positive, could be negative for large |A;|, 
implying instability. Without random motion, J- = 1 and 
short waves with k > fccrit or A < Acrit = 47r^G'I]/K^, will 
be unstable. Unlike for the Jeans instability in a stationary 
medium, longer waves are stabilized by rotation, embodied 
in the term. 

For disks with random motion, the reduction factor T — >■ 
for large irrespective of the frequency, for the reason 
given above. Thus random motion stabilize s short waves, 
as for the Jeans instability. iToomrd ()1964[ ) showed for a 
Gaussian distribution of radial velocities that the RHS of 
eq. (IS]) is > for all |fc| when an > (jR^ciit, where 

3.3601] , 

Cit.crit — • (11) 

K 

Thus Toomre's local ELxisymmetric stability crite- 
rion is Q = (Tjf/an^ciit > 1- When reasonably constant 
with radius, the locally estimated Q value is a good in- 
dicator of global axisymmetric stability. For exa mple, the 
local v alues in the equilibrium models proposed bv lKalnai^ 
are in reasonable agreement with those he derived 
from global axisymmetric modes (his Fig. 2). It should 
be noted that an axisymmetrically stable disk, i.e. for 
which Q > 1 everywhere, could still be unstable to non- 
axisymmetric modes. 

The local axisymmetric stability criterion for rotating 
gas disks is similar. The longest unstable wavelength, Acrit 
in a cold {vg = 0) disk is the same as for a stellar disk. 
The minimum of the quadratic expression in k on the RHS 
of eq. (|10l) occurs for k = TrGS/w^, and the condition of 
marginal stability (w^ = 0) at this minimum is readily 
solved, yielding t;.,.r.rit = ttGE/ k. 

Vandervoort ([19701 ) . iRomeol (|l992l ). and others have ex- 
tended the discussion of axisymmetric stability to disks of 
finite thickness. The principal difference is the dilution of 
the self-gravity term caused by spreading the disk matter 
in the vertical direction, resulting in a somewhat more sta- 
ble disk. The precise correction depends on the assumed 
vertical structure of both the disk and the wave, but is 
minor when the characteristic disk thickness zq <^ Acrit, 
which is usually the case. 

More realistic composite disks consist of multiple stel- 
lar populations having a range o f velocity dispersi o ns, as 
well as a cool gas component. iJog fc SolomonI (|l992f l 
considered the stability of two gravitation ally coup l ed gas 
disks having different sound speeds, but iRafikovl (|200lh 
improved their analysis to include multiple stellar compo- 
nents. These analyses show that while the gas component 
may contain a small fraction (typically < 10%) of the total 
mass, its low sound speed causes it to have a dispropor- 
tionate effect on the overall stability. 

3.4. Self- consistent Density Waves 

Fig. |4] gives, for four different values of Q, a graphical 
representation of the stellar (solid curves) p and gaseous 

^ The reduction factor used assumes a Gaussian distribution of radial 
velocities 
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Fig. 4. — The solid lines sho w the dispersion relation (|9)l, while 
the dashed lines are from l|10p . each for several values of Q. The 
ordinate s = \u — mf!|/K. Since the signs of both the abscissa 
and the ordinate are suppressed, the graph is reflection symmetric 
about both axes. Note that s = at CR and s = 1 at the Lindblad 
resonances. 

(dashed curves) dispersion relations eqs. ([9]) & (fTOl) re- 
spectively. The curves do not differ for Q = for which s 
is imaginary when fc > fccrit- The dashed and solid curves 
have similar forms for moderate fc but differ substantially 
for large fc, where frequencies in the stellar sheet cannot 
exceed |a; — m£l\ — k. Thus waves in a stellar sheet extend 
only from CR, where s = 0, to the Lindblad resonance on 
either side, where s = 10 

A gas layer, on the other hand, can support gravitation- 
ally modified sound waves of arbitrarily high frequency. 
Note that this important difference from the stellar case 
can lead to artificial modes when hydrodynamic equations 
are used as a simplifying approximation for a stellar disk; 
while meaningful results can be obtained using this ap- 
proximation, it is important to check that the derived 
modes do not cross Lindblad resonances. 

Since the dispersion relation (|9]) does not depend on the 
sign of k, curves for negative k, or leading waves, are sim- 
ple reflections about the k — Q axis. Furthermore, the 
relation gives the square of (w — mfi), implying that so- 
lutions when this frequency is negative, which arise inside 
CR, are simple reflections about the s = axis. Thus 
were the signs not suppressed. Fig. |4] shows only the panel 
for trailing waves outside CR, while the other quadrants 
would be reflected images. 

In the quadrant shown, there are either two values of k 
for each value of s, known as the short- and long- wave 
branches, or there are none. Only for the marginally 
stable case of Q = 1 do solutions exist for all frequencies 
< s < 1 and all |fc|. When Q > 1, a forbidden region 
opens up in the vicinity of CR where waves are evanescent. 

Recall that the relation (jH]) was derived assuming both 
\kR\ ^ 1 and that the waves are tightly wrapped ~ i.e. 
that the wave vector is directed radially. These approxima- 
tions may not be too bad on the short-wave branch (large 
k) but long- wave branch solutions may not exist at all, ex- 
cept in a very l ow-mass disks where fccri t is lar ge, and must 
fad as ^ 0. iLvnden-Bell fc KalnaisI (119721 ) gave an im- 

Additional Bernstein- type waves exist near integer values of s > 1 , 
but such solutions seem to be of little dynamical importance. 



proved dispersion relation (their eq. All) that relaxed the 
requirement that the wave vector be radial, but still used 
the plane wave potential ([7]) that requires |fci?| ^ 1. Their 
relation is no different for tightly-wrapped waves, but has 
no long-wave branch or forbidden region for open waves 
when Q > 1. 

3.5. Group Velocity 

The waves in a disk described by the local dispersion 
relations ^ & (fTOj) have a phase speed, equal to the pat- 
tern s peed, in the azimuthal direction. However, iToomrj 
(|1969f ) pointed out that these dispersive waves also have a 
group velocity that is directed radially. Since Vg — duj/dk, 
the group velocity is proportional to the slope of the lines 
in Fig. |4l Each portion of a wavetrain, or wave packet, is 
carried radially at the group velocity, retaining its pattern 
speed u whilst gradually adjusting its wavelength, and also 
transporting wave action or angular momentum. 

For trailing waves outside CR, the situation illustrated 
in Fig. HJ the group velocity is positive on the short-wave 
branch, and the waves carry angular momentum outwards 
towards the OLR. The curves have the opposite slope in- 
side CR, where density waves are disturbances of nega- 
tive angular momentum (s tars there lose angular momen - 
tum as the wave is set up, iLvnden-Bell fc Kalnaisl[l972D . 
and therefore the inward group velocity on the short-wave 
branch leads once more to outward angular momentum 
transport. Outward angular momentum transport is in the 
sense e xpected from the gravitationa l stresses of a trailing 
spiral (|Lvnden-BeU fc Kalnaislll972f ). 

However, the sign of the slope du/dk on the dubious 
long-wave branch is opposite to that of short waves, in- 
dicating that angular momentum in that regime is trans- 
ported in the direction opposite to the gravity torque! This 
apparent contradiction is resolved by the phenomenon of 
lorry transport, a term coined bv lLvnden-Bell fc Kalnai j 
(|l972il to describe an advective Reynolds stress (see also 
Appendix J of BT08). Stars in their epicyclic oscillations 
gain angular momentum from the wave near apo-center 
and lose it back to the wave near peri-center; no star suf- 
fers a net change, but angular momentum is still carried 
inward at a sufficient rate to overwhelm the gravity torque. 

Thus where eq. Q holds in a disk with approximately 
uniform Q > 1, a tightly wrapped trailing wave packet 
originating on the long-wave branch, will travel towards 
CR until it reaches the edge of the forbidden zone where 
it "refracts" into a short wave that carries it back towards 
the Lindblad resonance. The details of the turning point 
at the forbidden zone, which requires an evanescent wave 
propagating into the forbi dden region, were descr ibed by 
Mark (,19_76,) and bv Goldreich fc Tremaii^ (|l978|). 

The fate of the short-wave propagating towards the 
Lindblad resonance also requires more sophisticated anal- 
ysis because the dispersio n relation © does not hold near 
resonances. iMarkI (jl974f ) carried the analysis to second 
order to show that the wave is absorbed there through 
the wa ve-particle coupling described bv Lvnden-Bell & 
Kalnajsll972D. 

3.6. Swing Amplification 

Goldreich fc Lvnden-BeU (|1965b[). rJulian fc Tooinri (|T966h 
and iToomrg (,1981.) extended the above analysis for tightly- 
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wrapped waves to waves of arbitrary inclination. They 
again adopted the approximate potential of a short- wave- 
length plane wave (eq. [7]), and focused on a local patch of 
the disk whose dimensions were assumed small compared 
with Rq, the radial distance of the patch center, which or- 
bits the disk ceii t er at t he local circular speed V — Rq^Iq. 

Following iHilll (jl878f ). they then introduced Cartesian- 
like coordinates, x = R — Rq and y = Ro{4> — 4>q), where 
(pQ is the azimuth of the orbiting center of the patch. The 
vertical coordinate, z, is unchanged from the usual cylin- 
drical polar coordinate. Since the angular rotation rate 
decreases outwards in galaxies, the motion of stars in the 
patch is that of a shear flow, with a speed that varies as 
y = —2Ax, together with Coriolis forces arising from the 
rotation of the patch at angular rate JIq. The Oort "con- 
stant" A = —^Rdfl/dR. The undisturbed surface density 
of matter, S, is assumed constant over the patch. This set 
of approximations is described as the sheared sheet (see 
BT08 pp. 678-681). 

If a is the inclination angle of a plane density wave to 
the azimuthal direction, with 180° > a > 90° for leading 
waves and 90° > a > 0° for trailing waves, then the shear 
causes the pitch angle to change with time as cot a — 2At, 
with t = when a = 90°, i.e. when t he wav e is purely ra- 
dial. Bot hlGoldreich fc Lvnden-Belll (I1965H ) for a gaseous 
sheet and'Julia n fc Toomrd (jl966'^ for a stellar sheet found 
waves amplify strongly, but transiently, as they shear from 
leading to trailing. 

The unremitting shear flow ultim ately tears apart any 
transient disturbance, and indeed Jul ian &: Toomrd (jl966f ) 
asserted that the stellar sheet is locally stable at small- 
amplitudes whe n Q > 1, but th ey did not give the proof. 

Fig. \E\ from iToomrd ()1981h . illustrates the dramatic 
transient trailing spiral that results from a small input 
leading disturbance. This calculation illustrates that even 
a global, small-amplitude disturbance does not persist, but 
decays after its transient flourish. In the late stages, the 
"wave action" propagates at the group velocity ( ^3.5p away 
from corotation, where it is absorbed at the Lindblad res- 
onances. 

Notice also that the wave in this Q = 1.5 disk extends 
across the forbidden region near CR. The absence of solu- 
tions near s = when Q > 1 in Fig. |4]is a consequence of 
assuming a steady, tightly-wrapped wave; shearing waves 
do not respect this restriction. 

The responsiveness of a disk to input leading signal de- 
pends both on the degree of random motion - it decreases 
rapidly over the range 1 < Q < 2 - and upon the az- 
imuthal wavelength, which i s A,, = 2tt R/m for an m-armed 
pattern. lJulian fc Toomrd (|1966l ) defined the parameter 



X 



Rk 



crit 



Acrit m 

In a disk having a flat rotation curve, iToomrd (|l98ll ) 
showed that the strongest amplification occurs when 1 < 
X < 2.5, or about twice the scale of the longest unsta- 
ble axisymmetric wave. The peak amplification shifts to 
larger X values for declining rotation curves and vice versa 
for rising ones. Naturally, the phenomenon of swing am- 
plification disappears entirely in uniformly rotating disks. 

By making the assumption that spiral patterns in galax- 
ies are amplifie d solely bv this mecha nism. Athanassoula 
et al~ Jl987^ argued that the multiplicity of the observed 







Fig. 5. — The time evolution of an input leading wave packet in 
the half-mass Mestel disk. The unit of time is half a circul a r orb it 
period at the radius marked CR. Reproduced from I Toomrd l|198ir . 



spiral arm pattern can be used to place bounds on the mass 
density in the disk. The circular orbit speed in a nearly ax- 
isymmetric galaxy affords a direct estimate of the central 
attraction, but it is hard to determine how much of the to- 
tal attraction is contributed by the disk, and how much by 
other components, such as the dark matter halo. Because 
Acrit is proportional to the disk surface density, for a given 
rotation curve, the most vigorous amplification shifts to 
higher sectoral harmonics in lower mass disks. Generally, 
disks that contribute most of the central attraction, i.e. 
close to mELximum disks, would prefer bi-symmetric spi- 
rals, while higher arm multiplicities would indicate signif- 
icant ly sub-maximal disks (see also ISellwood fc Carlbergj 
[1981 . 

4. BAR INSTABILITY 

The puzzle of how galaxy disks could be stable presented 
a major obstacle to the development of our understanding 
of disk dynamics for many years. Superficially reasonable 
models of dis k galaxies wer e found repeatedl y, both in sim- 



ulatio ns fe.q.lMiller et al\ K970: Hohl 19 7ll: IZang fc Hohll 



119781: ICombes k SanderdlT981.: ,Sellwoo(Jll98^ Athanas- 



soul a fc Sellwood 119861: iKhoperskov et aLH2007l: Dubinski 
et al. 12009^ and in linear stability analyses (KalnajiHazi; 
lAoki e t aj, 1979; Toomrc 1981; Sawamura,19 88: Vauterin 
fc Deionghe 119961: iPichon fc CannonI Il997l: Korchagin et 
al. 120051 : lPolvachenkoll2005l : lJalalill2007f r to possess vigor- 
ous, global and disruptive bar-forming instabilities. Fig. [5] 
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Fig. 6. — The formation of a bar in an (unpublished) Af-body 
simulation. The three orthogonal projections show only particles in 
the initially exponential disk, halo particles are omitted. The disk, 
which started with Q = 1.2, has unit mass and unit scale length and 
G = 1. One orbit at ij = 2 takes 15 time units, where the central 
attractions of the disk and halo are nearly equal. 



illustrates the global, disruptive nature of this instability. 
While it is premature to claim that this problem has been 
completely solved, it now seems that the stability of a 
massive disk galaxy requires only a dense bulge-like mass 
component near the center, and owes little to the inner 
density of a dark matter halo. Galaxies lacking a central 
mass concentration, however, are still believed to require 
significant inner halo mass for global stability. 

4.1. Mechanism for the Bar Mode 

Toomre (jl98ll ) provided the most important step for- 
ward by elucidating the mechanism of the bar instability 
(see also BT08 §6.3). Linear bar-forming modes are stand- 
ing waves in a cavity, akin to the familiar modes of organ 
pipes and guitar strings. Refiections in galaxies take place 
at the center and at the corotation radius, except that out- 
going leading spiral waves incident on the corotation cir- 
cle are super-reflected into amplified ingoing trailing waves 
(i.e. swing amplification), while also exciting an outgoing 
transmitted trailing wave. The feedback loop is closed by 
the ingoing trailing wave reflecting off the disk center into 
a leading wave, which propagates outwards because the 
group velocity of leading waves has the opposite sign to 
that of the corresponding trailing wave. The amplitude of 
the continuous wave-train at any point in the loop rises 
exponentially, because the circuit includes positive feed- 
back. 

Toomre supported this explanation with linear stabil- 
ity studies of two disks. The Gaussian disk, which has a 
low central density, manifests a set of modes in which the 
more slowly- growing, higher "overtones" display the kind 
of standing wave pattern to be expected from the super- 
position of ingoing trailing and outgoing leading waves. 



(The eigenfrequencies of these modes are shown below in 
Fig.l) 

The other linear stability study he presented was for the 
inappropriately-named "Mestel" diskj^ whose unusual sta- 
bility prope rties were flrst described by Toomre's student 
IZand (il976D in. an unpublished thesis, and later bv Evans 
& Read (|1998( ). This disk has the scale-free surface den- 
sity profile E = /{2ttGR), with Vc being the circular 
orbital speed that is independent of R. Zang had carried 
through a global, linear stability analysis of this disk, with 
random motions given by a smooth distribution function. 
In order to break the scale- free nature of the disk, Zang in- 
troduced a central cutout, and later an outer taper, in the 
active mass density, replacing the removed mass by rigid 
components in order that the central attraction remained 
unchanged at all radii. The dominant linear instabilities 
he derived for the tapered disk were con firmed in A^-body 
simulat i ons bylSellwood fc EvansI (|200l[ ). 

Zang (jl976i ) showed that a full-mass Mestel disk is stable 
to bi-symmetric mo des, even if Q < 1 , prov ided the tapers 
are gentle enough. lEvans fc ReadI (|1998( ) extended this 
important result to other power-law disks, finding they 
have similar properties. Such disks are not globally stable, 
however, because they suffer from lop-sided instabilities 
(see ^5^. By halving the active disk mass, with rigid mass 
preserving the overall central attraction, iToomrd (|l98lh 
was well able to eliminate the lop-sided mode from the 
Mestel disk. In fact, he claimed that when Q = 1.5, the 
half-mass Mestel disk was globally stable, making it the 
only known model of a non-uniformly rotating disk that is 
stable to all small amplitude perturbations. 

Cutting the feedback loop really does stabilize a disk 
(see next section), which seems to co nfirm Toomre ' s cav- 
ity mode mechanism. Despite this, iPolvachenkol (|2004l ) 
pointed out that the strong emphasis on phenomena at 
corotation places the blame for the instability on a radius 
well outside where the resulting bar has its peak ampli- 
tude. He therefore proposed an alternative mechanism for 
the bar instabilitv based u pon orbit alignment f Lvnden- 
Bell"Ii979, see also TO.ip . Even though the mechanism 
was ori ginallv envisaged as a slow trapping process. Polv- 
achenko (|2004l ) argued it may also op erate on a dynamical 
timcscale, and he devised (jPolvachen ko 2005) an approx- 
imate technique to compute global modes that embodied 
this idea. While his method should yield the same mode 
spectrum as other techniques, his alternative characteriza- 
tion of the eigenvalue problem may shed further light on 
the bar-forming mechanism. 

4.2. Predicted Stability 

The fact that small-amplitude, bi-symmetric instabili- 
ties are so easily avoided in the Mestel disk, together with 
his und erstanding of bar-forming modes in other mod- 
els, led IToomrd ()198lD to propose that stability to bar- 
formation merely required the feedback cycle through the 
center to be cut. He clearly hoped that a dense bulge- like 
mass component, which would cause an ILR to exist for 
most reasonable pattern speeds, might alone be enough to 
stabilize a cool, massive disk. 

Unfortunately, this prediction appeared to be contra- 



^ IMestell l l 19631) solved the far greater challenge of a disk of finite 
radius that has an exactly flat rotation curve. 
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di eted almost immediately by the findings of Efstathiou et 
al. (| 19821 ) ■ whose -/V-bodv simulations formed similar bars 
on short timescales irrespective of the density of the central 
bulge co mponent! They seemed to confirm previous con- 
clusions (jOstriker fc Peebleslll973f ) that only significantly 
sub-maximal disks can avoid disruptive bar-forming insta- 
bilities. 

However. ISellwoodI (llQSQa*) found that Toomre's predic- 
tion does not apply to noisy simulations because of non- 
linear orbit trapping. Simulations in which shot noise was 
suppressed by quiet start techniques did indeed manifest 
the tendency towards global stability as the bulge was 
made more dense, as Toomre's linear theory predicted. 
Shot noise from randomly distributed particles, on the 
other hand, produced a large enough amplitude collective 
response for non-linearities to be important, and a noisy 
simulation of a linearly stable disk quickly formed a strong 
b ar consistent with the results reported bv Efstathiou et 
al. (|1982[ ). 

Since density variations in the distribution of randomly 
distributed particles are responsible for bar formation in 
this regime, the reduced shot noise level from larger num- 
bers of particles must result in lower amplitude responses 
that ultimately should avoid non-linear trapping. The pre- 
cise particle number required depends on the responsive- 
ness of the disk, which is weakened by random motion, 
lower surface d ensity, by increased disk thickness, or grav- 
ity softening. lEfstathiou et al\ (|1982f ) employed merely 
20 000 particles, which was clearly inadequate to capture 
the linear behavior. Howeve r, robustly stable, massive 
disk s have been simulated b y ISellwood fc Moor3 (|l999f ) 
and lSellwood k, EvansI (|200l[ ) that employed only slightly 
larger particle numbers. Note that these latter models also 
benefitted from more careful set-up procedures to create 
the initial equilibrium. 

Thu s the stabilizing mechanism proposed by iToomrd 
(|l98lD does indeed work in simulations of high enough 
quality, and presumab l y ther efore also in real galaxy disks. 
Indeed, iB arazza et al\ (120081 ) found a decreased incidence 
of bars in galaxies having luminous bulges, and argued that 
their result supports Toomre's stabilizing mechanism. 

Thus the absence of bars in a significant fraction of high- 
mass disk galaxies does not imply that the disk is sub- 
maximal. The old stability criteria proposed by Ostriker & 
Peeb les (1973). Efstathiou et al. (.19821 and Christodoulou 
et al. (1 19951) apply onl y to d isks that lack dense centers; 
indeed lEvans fc ReadI (|1998[ ) showed explicitly that the 
power-law disks are clear counter-examples to the simple 
bi-svm metric stability criterion proposed by Ostriker & 
PeeblesHiza!). 



4.3. Residual Concerns 

While all this represents real progress, a few puzzles re- 
main. The most insistent is the absence of large, strong 
bars in galaxies like M33, which has a gently rising rota- 
tion curve. A hhough many spiral arms can be counted in 
blue images (ISandage fc Humphreysl [T980[ ) . the near IR 
view ([Block et aLll2004l ) manifests an open two-arm spiral, 
suggesting that the disk cannot be far from maximal, and 
also reveals a mil d bar near the center of the galaxy. Cor- 
belli & Walterbos ()2008f ) also found kinematic evidence 
for a weak bar. Perhaps the stability of this galaxy can 



be explained by some slightly larger halo fraction, or per- 
haps the weak bar has some unexpected effect, but there is 
no convincing study to demonstrate that this galaxy, and 
others like it, can support a 2-arm spiral without being 
disruptively unstable. 

The second concern is that lop-sided instabilities appear 
in extended full-mass disks with fiat or declining rotation 
curves, which is discussed next. A third concern, which 
is discussed in !j9.2[ is that the mechanism is unable to 
predict the presence or absence of a bar in a real galaxy. 

5. LOP-SIDED MODES 

Many galaxies have apparently lop-sid ed disks. The 
treatment here will not go into detail, since I Jog fc CombesI 
(120091 ) have recently reviewed both the observational data 
and theory. 

Both theoretical and simulation work on m = 1 distor- 
tions to an axisymmetric disk require special care, since 
the absence of rotational symmetry can lead to artifacts 
unless special attention is paid to linear momentum con- 
servation. Rigid mass components present particular diffi- 
culties, since they should not be held fixed, and extensive 
mass components are un li kely t o respond as rigid objects. 

As noted above, IZangj (jl976[ ) found that the dominant 
instability of the centrally cut-out Mestel disk was not the 
usual bar instability, but a lop-sided mode, which persists 
in a full mass disk no matter how large a degree of ran- 
dom motion or gentle the cutouts. This surprising finding 
was confirmed a nd ex tended to general power-law disks by 
I Evans fc R ead (1998'). A lop-sided instability dominated 
simulations (ScUwood 1985) of a model having some re- 
semblance to Zang's, in that it had a den s e mas sive bulge 
and no extended halo, while ISaha et all (|2007D reported 
similar behavior in sim ulations of a bare exponential disk. 
iLovelace et al\ ()1999l ) found pervasive lop-sided instabili- 
ties near the disk center in a study of the collective modes 
of a set of mass rings. 

Various mechanisms have b een proposed to account for 
this instability . Baldwi n et al. ([l980) and Earn fc Lynden- 
Bell (|1996^ explored the idea that long-lived lop-sidedness 
could be constructed from cooperative orbital responses of 
the disk stars, along the lines discussed for bars bv Lynden- 
BeU (I1979D . iTremaind (I2005D discussed a self-gravitating 
secular instability in near-Keplerian potentials^ A more 
promising mechanism is a c avity mode, simil ar to that for 
the bar- forming instability (|Dury et aLll2008D : the mecha- 
nism again supposes feedback to the swing- amplifier, 
which is still vigorous for m = 1 in a full-mass disk. 

Feedback through the center cannot be prevented by an 
ILR for m = 1 waves, since the resonance condition f2p = 
— (eq.[S]) can be satisfied only for retrograde waves. 
But the lopsided mode can be stabilized by reducing the 
disk mass, which reduces the X par ameter (eq. [T^ until 
amplific ation dies for m = 1 (iToom rc 1981). ScUwood 
& Evans (|200ll ) showed that, together with a moderate 
bulge, the dark matter required for a globally stable disk 
need not be much more than a constant density core to 
the minimum halo needed for a flat outer rotation curve. 

A qualitatively different lop-sided instability is driven 



^ The "sling amplification" mechanism proposed bv lShu et al\ 1119901) 
applies only to gaseous accretion disks, since it relies on sound waves 
propagating outside the OLR. 
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Fig. 7. — The top panel illustrates a groove in the sheared sheet 
f i|3.6l l model. The light shaded region has the full undisturbed disk 
surface density, S, which in this frame has a linearly varying shear 
flow. The white region is the groove, which has lower surface density, 
and the density excesses brought into the grove by small wave-like 
disturbances on either side of the groove are highlighted by the dark 
shading. The width of the groove is exaggerated for clarity. 
The lower panel, which has the correct aspect ratio, shows the sup- 
porting responses to the density variations on the groove edges. The 
dashed lines mark the locations of the Lindblad resonances. 



by counter-rotation. Th is second kind of m = 1 instabil- 
ity was first reported by Zang fc Hohl (|l978f ) in a series of 
A^-body simulations designed to explore the suppression of 
the bar instability by reversing the angular momenta of a 
fraction of the stars; they found that a lop-sided instability 
was aggravated as more retrograde stars were included in 
their attempts to subdue the bar mod e. Analyses of vari- 
ous disk inodels with retrograde stars (fArakilll987l: Sawa- 
mura ]l988l : lDurv et aLll2008l ) have revealed that the growth 
rates of lop-sided instabili ties increase as the fracti on of 
retrograde stars inc r eases . iMerritt &: Stiavellil ()199Clf ) and 
ISellwood fc Valluril (|1997f ) found lop-sided instabilities in 
simulations of oblate spheroids with no net rotation. Their 
flatter models had velocity ell ipsoids with a strong tangen- 
tial bias, whereas ISellwood fc Mcrritt (1994) found that 
disks with half the stars retrograde, together with moder- 
ate radial motion w ere surprisingly stable. 

Weinberg ()1994f ) pointed out that lop-sided distortions 
to near spherical systems can decay very slowly, leading 
to a protracted period of "sloshing". This seiche mode 
in a halo, which decays particularly slowly in mildly con- 
centrated spherical s ystems, could pr ovoke lop-sidedness 
m an embedded disk (jKornreich et al. 2002; Idcta 2002). 

If lop-sidedness is due to instability, then the limited 
work so far suggests that it would imply a near-maximum 
disk. But lop-sidedness in the outer parts could also be 
caused by tidal interactions, or simply by asymmetric disk 
growth, with the effects of differential rotation being mit- 
igated pe rhaps by the cooperative o rbital responses dis- 
cussed bv lEarn fc Lvnden-Belll (|1996[ ). 

6. GROOVE AND EDGE MODES 



Not all true global instabilities in disks require a feed- 
back loop. Another class of modes is simply driven from 
corotation, with trailing waves propagating both inwards 
and outwards to be absorbed, in stellar disks, at the Lind- 
blad resonances on either sid e. 

Lovelace & Hohlfeld (|l978f ) showed that disks are desta- 
bilized by a local extremum in the radial variation of the 
ratio of surface density to vorticity, S/|Vxu| = I]ri/(2K^), 
or the reciprocal of potential vorticity, and proposed that 
the instability created flat rota tion curves for wh i ch the 
potential vorticity is also flat. ISellwood &: Kah"nl (|1991h 
demonstrated that the instability caused by a rather 
insignificant, narrow, decrease in surface density, i.e. a 
"groove" , is a global spiral mode. 

The mechanism is easiest to visualize in a disk without 
random motion, where small-scale surface density varia- 
tions are not blurred by epicyclic motions. In this case, a 
deficiency of stars over a small range of angular momen- 
tum creates a groove in the surface density, as shown in the 
sheared sheet model in the top panel of Fig. [71 The groove 
itself is unstable because of the gravitational coupling be- 
tween disturbances on either side. The dark shaded areas 
in the Figure illustrate regions where small sinusoidal ra- 
dial displacements of material on each edge have created 
high density regions where the density was previously low. 
Disturbing gravitational forces arise from the density ex- 
cesses, as illustrated, which are directed along the groove 
if the wavelength is long compared with the groove width. 
Material that is pulled back loses angular momentum and 
sinks toward the center of the galaxy, while that which is 
urged forward gains and rises outwards. Thus each den- 
sity excess pulls on the other across the groove in such a 
way as to cause it to grow exponentially, i.e. the combined 
disturbance on the two sides is unstable. The groove edges 
need be only steep gradients, not discontinuities. Further- 
more, the mechanism is the same, but harder to visualize, 
in a disk with random motion where the density of stars 
is depleted over a narrow range of angular momentum. 

The growing disturbance in the groove creates wave- 
like mass variations along the groove that are effectively 
growing co-orbiting mas s clumps of the type envisaged by 
lJulian fc Toomrd (|l966f ). Generalizing their apparatus to 
allow for exponentially growing masses. ISellwood &: KahnI 
(1991') estimated the expected disk response, as shown in 
the lower panel of Fig. [Tjfor the parameters Q = 1.8 and 
X — 2, and low growth rate. The disk supporting response 
transforms the quite trivial disturbance in the groove into 
an extensive spiral instability! 

Unlike the sheared sheet, the azimuthal wavelength in 
a full disk can take on only discrete values, Xy = 
which are all unstable, but the one which grows most 
rapidly is that for which swing-amplification (i j3.6p is the 
most vigorous, i.e. 1 < AT < 2.5 (eq. [T^ in a flat rotation 
curve. The simulation shown in Fig. |8] illustrates the scale 
and vigor of the mode, which saturates at ^ 20% overden- 
sity because of the onset o f horseshoe orbits at corotation 
(Sellwood & Binnev 2002). 

Lovelace & Hohlfeld (197j) and lSellwood fc KahnI (fT99l[ ) 
found that almost any narrow feature in the angular mo- 
mentum density is destabilizing, although the modes of a 
simple ridge, for example, come in pairs with CR some 
distance from the ridge center. 
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Fig. 8. — The later part of the growth and subsequent decay of an isolated s piral mode in a disc that w as seeded with a groove. Disk 
rotation is counter-clockwise and disturbance forces in this simulation, taken from lSellwood fc BinnevI II2002I) . were restricted to m = 2 only. 




Fig. 9. — The variation of pattern speeds and growth rates for 
the first few m = 2 mo des of the Gaus sian disk as the active disk 
mass is decreased. From iToomrel l|198ir) . 

Edge modes (|Toomrelll981l : iPapaloizou fc Liij[T989h 

are close cousins of groove modes, and the mechanism can 
be understood in the same fashion. Density variations 
on a single edge without a supporting disk response are 
neutrally stable, but the necessarily one-sided wakes in 
the massive interior disk add angular momentum to the 
density enhancements on the edge, causing them to grow. 

The eigenfrequencies of the low order bisymmetric 
modes of the Gaussian disk vary w ith decre a sing d isk mass 
as shown in Fig. O taken from iToomrd (jl98l[ ). Most 
are cavity modes, but the edge-mode, labelled "mode D" 
stands out because its growth rate in particular declines 
more slowly with decreasing disk mass. If the mechanism 
for both types of mode involves swing amplification, one 
might expect their vigor to be similarly affected by the rise 
in the X parameter as the surface density decreases. But 
the decreased surface density also slows the group velocity 
while the rapidly declining pattern speed of each cavity 
mode moves CR farther out in the disk. Thus the growth 
rates of cavity modes drop more quickly because of the 
increased travel time for the wave packet to complete the 



feed-back loop, a factor that is absent for the edge mode. 

The edge instability requires only a steep gra dient in 
the surface density, which need not drop to zero. iToomrel 
(fT989il gave the condition for instability as "the radial dis- 
tance over which the disc density undergoes most of its 
rapid change should be no larger than about one quar- 
ter of the axisymmetric stability length Acrit" , assuming 
"the disc is massive and cool enough for vigorous swing- 
amplification." 

Curiously, global bar-forming modes were originally 
thought to be related the rotational instabilities of uniform- 
d cnsity Maclaurin sphe roids of incompressible fluid 
(jOstriker fc Peebleill973f ) , an idea reinforced by the vigor 
of th e bisymmetric mode of the sha rp-edged Maclaurin 
disk (|Kalnaisl[T97l . IToomrel (119811 ) noted that the in- 
stability of such unrealistic galaxy models may be more 
closely related to the edge mode than to the cavity mode 
described in 

As the consequence of an edge instability in a realistic 
galaxy model is to blur the edge, it seems unlikely that 
galaxy disks can retain unstable density gradients for in- 
terestingly long periods. Similarly, the large-amplitude 
evolution of groove or ridge modes quickly erases the fea- 
ture that gave rise to them. Nevertheless, the modes have 
other possible consequences (see ^7.6^ . 

7. SPIRAL STRUCTURE 

Despite many decades of effort, no theory to account 
for the graceful spiral patterns in disk galaxies is widely 
accepted. Most workers in this field agree that spiral pat- 
terns are gravitationally driven variations in surface den- 
sity in the old stellar disk, as s upported by photom etric 
data {e.g. lGrosb0l "all l2004l : IZibetti atl I2009D and 
strea ming motions in hi gh spatial resolution velocity maps 
{e.g. IShettv et ai]\2007\ ). 

There seems litt le doubt that some spiral patterns are 
tidally driven fe.o. ISalo fc Laurikainen|[l993l : iDobbs et al\ 
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|2010[). while othe rs could be the driven responses to bars 
( Buta et aI]|2009D . Although these two ideas may account 
for a large fraction of the cases (jKormendv fc NormarJ 
I1979I ) , es pecially if orbiting da rk matter clumps can excite 
patterns (jPubinski et a/.l[2008l ) , spirals can still develop in 
the absence of either trigger, as revealed in simulations. 

The idea that spirals could be self-excited oscillations 
of the stellar disk represents the greatest theoretical chal- 
lenge. While there is general agreement that gas seems to 
be essential, no picture seems complete, and current the- 
ories disa gree on even t he lifetimes of the patterns. One 
idea (e.g. iBertin fc Linl 1996), is that spiral features are 
manifestations of quasi-steady global modes of the under- 
lying disk. Alternatively, they could be short-lived, recur- 
rent transient patterns that originate eithe r from forcing 
by density fluctuations {e.g. Toomre '1990"), or else from 
recurrent vigorous instabilities {e.g. ScUwood 2000). 

A serious obstacle to progress in this area has been the 
absence of observational discriminants that would favor 
one of these radically differing viewpoints over the other. 
The predictions for density variations or gas responses at a 
single instant are essentially independent of the generating 
mechanism a nd do not dep e nd st rongly on the lifetime of 
the pattern. iMeidt et al\ (|2009l ) attempted to measure 
radial variations of patter n speeds using a generalization 
of the method devised by iTremaine Sz WeinbergI ()l984al ) 
(see t|9.5p . They reported lower pattern speeds at larger 
radii, but their measurements still tell us little about spiral 
lifetimes. However, the velocity distribution of stars in the 
solar neighborhood ()Nordstrom et a/.ll2004l . and Fig. [T]) is 
most naturally accounted for in the transient picture (see 

7.1. Spirals as Global Modes of Smooth Disks 

Simple mode ls of disk ga laxies possess many linear in- 
stabilities {e.g. lJalalill2007f ). The bar-forming mode O 
is usually the fastest growing, but it has almost no spiral- 
ity. These studies are therefore important to understand 
stability, but even the higher modes are not particularly 
promising for spiral generation^ 

The "den sity wave" theory fo r spiral modes, described 
in detail by iBertin fc Linl ()1996( ) , invokes a more specific 
galaxy model with a sub-maximal disk that is dynamically 
cool in the outer parts and hot in the inner disk. The 
local stability parameter is postulated to be 1.0 < Q < 
1.2 over most of the disk and to rise steeply to Q > 2 
near the center. Bertin, Lin, and their co-workers perform 
a global analysis using the hydrodynamic approximation 
(BT08 §5.1), which reveals slowly growing spiral modes 
under these specific conditions. 

The mechanism (Mark 1977) is a cavity mode, having 
qualitative similarities to that for the bar mode (S3Ji 
but the tightly- wrapped waves are trailing around the en- 
tire cycle. The inner turning point is at a Q-barrier: a 
steeply increasing Q value causes the forbidden zone (§ 



^ IKorchagin et oil ([2009) calculated essentially gas-dynamical modes 
for models of specific galaxies, and argued that the shapes of one, 
or more, of the lower-order modes could be matched to the observed 
spiral pattern. Ifowever, it is unclear that rapidly growing, linear 
modes can be seen for long at finite amplitude before non-linear 
effects will change their appearance, and it seems even less likely 
that two modes with different growth rates should have similar large 
amplitudes at the time a galaxy is observed. 



to broaden to the point that in-going short waves get "re- 
fracted" into outgoing long waves, which prevents the wave 
train from reaching the ILR where it would be damped. 
The long waves then propagate out to near the CR in 
the Q ^ 1 part of the disk, where they switch back to 
the shortwave branch wit h a small de gree of amplification. 
The WASER mechanism ()Marklll976h at this turning point 
involves a third, transmitted wave that is "radiated" out- 
wards on the short-wave branch, carrying away the angular 
momentum to excite the mode in the inner disk. Thus the 
amplifier involves a small "swing" from the long- to the 
short-wave branch, whereas the bar instability uses a full 
swing from leading to trailing. This difference, together 
with their assumption that the disk is sub- maximal {i.e. 
X > 3 fo r TO = 2, see eg. [T^. allows the mode to be slowly 
growing. iLowe et al\ ([l994) present a model of this kind 
to account for the spiral structure of M81. 

In order to justify the "basic state" of the disk they re- 
quire, [Bertin^Li^ ([iMl) argued heuristically that rapidly 
evolving features would have disappeared long ago and 
that low-growth-rate instabilities in a cool disk, created 
by gas dissipative processes and star formation, will dom- 
inate at later times. They invoked shocks in the gas to 
limit the amplitude of the slowly growing mode, leading 
to a quasi-steady global spiral pattern. 

The main objection to their picture is that it is likely 
that such a lively outer disk will suffer from other, more 
vigorous, collective responses with to > 2 that will quickly 
heat the disk, as described in H7.51 and destroy their pos- 
tulated background state. 

7.2. Recurrent Transients 



From the early work bvlMiller et al\ (IT970l') . iHohll 
iHockney fc Brownrigg (Il974f). James fc Sellwood 



1971 ) 



mm 



and Sellwood fc Carlberd ( 1984 ). iV-body simulations of 
cool, sub-maximal disks have exhibited recurrent transient 
spiral activity. This basic result has not changed for sev- 
eral decades as numerical quality has i mproved. 

Cla ims of long-lived spiral waves (e.g.YTh omasson et ali 
Il990f ) have mostly b een based on simulations of shor t du- 
ration. For example, lElmegreen fc ThomassonI (|1993D pre- 
sented a simulation that displayed spiral patterns for ~ 10 
rotations, but the existence of some underlying long-lived 
wave is unclear b ecause the pattern changed f rom snap- 
shot to snapshot. iDonner fc ThomassonI (jl994[ ) ran their 
simulations for fewer than 7 disk rotations and the bi- 
symmetric spiral they observed appeared to be an incip- 
ient bar instability. Zhang (1996,1998) adopted a similar 
mass distribution and also reported long-lived patterns in 
her simulations. The author has attempted to reproduce 
her results, and indeed obtained similar bi-symmetric fea- 
tures, but they appear to be the super-position of several 
waves having differing pattern speeds. 

Sellwood & Carlberg ()1984f ) stressed that patterns fade 
in simple simulations that do not include the effects of gas 
dissipation; the reason is the disk becomes less responsive 
as random motion rises due to particle scattering by the 
spiral activity f ^ilOl) . which is therefore self-limiting. They 
also demonstrated that mimicking the effects of dissipative 
infall of gas, such as by adding fresh particles on circular 
orbi ts, allowed patterns to recur "indefinitely." Later work 
(e.q. lCarlberg fc Freedmanlll985l : l^omrelll990D has shown 
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that almost any method of dynamical cooling can maintain 
spiral activity, as also happens in mo dern galaxy formation 
simulations fe.o. lAgertz et a/.ll201Cll ). 

Thus the transient spiral picture offers a natural expla- 
nation for the absence of spiral patterns in SO disk galaxies 
that have little or no gas; maintenance of spiral activity 
requires a constant supply of new stars on near-circular 
orbits. Other pieces of indirect evidence that also favor 
the transient spiral picture are reviewed in ijlOl 

7.3. Spirals as Responses to Density Fluctuations 

Goldreich & Lynden-Bell (ll965bD and iToomrd (|1990[ ) 

suggested that a large part of the spiral activity observed 
in disk galaxies is the collective response of the disk to 
clumps in the density distribution. As a spiral wake is the 
collective r esponse of a disk to an individual co-orbiting 
perturber (jJulian fc Toomr3 Il966f ). multiple perturbers 
will create multiple responses that all orbit at different 
rates. The behavior of this polarized disk reveals a chang- 
ing pattern of trailing spirals, which can equivalently be 
regarded as swing-am p lified (j j3.6|) noise^l 

I'oomre & Kalnajs ()1991h studied the amplified noise 
that arose in their A^-body simulations of the sheared sheet. 
The massive particles themselves provoke spiral responses 
with an amplitude proportional to the input level of shot 
noise, caused by density variations in the distribution of 
randomly distributed particles. Comparison of their ex- 
pectations with linear theory predictions revealed that the 
simulations were livelier than they expected, by a factor 
< 2, apparently from a gradual build-up of correlations 
between the mean orbital radii of the particles. 

This could be a mechanism for chaotic spirals in very 
gas rich discs, where a high rate of dissipati on may be able 
to maintain the responsiveness of the disc ()Toomrelll990l ) 
while the dumpiness of the gas distribution may make the 
seed noise amplitude unusually h igh. However, it se ems 
likely that spiral amplitudes (e.g. IZibetti et al\\200% are 
too large to be accounted for by this mechanism in most 
galaxies. Also, the spiral structure should be chaotic, with 
little in the way of clear symmetry expected. 

7.4. Non-linear Spiral Dynamics 

Ta gger and his co-workers (iTagger et a/.lll987l: Svgnet 
et g/. ll987l : lMasset fc Taggeijll997f ) suggested that global 
modes in stellar disks could be coupled through non-linear 
interactions. They proposed that wave 1 excites wave 2 
through second-order coupling terms that are large when 
CR of wave 1 lies at approximately the same radius as the 
ILR of wave 2. Conservation rules require a third wave 
such that 7713 = mi ±TO2; most likely an axisymmetric 
wave (ma = 0) if mi — m2. Many examples of multiple 
waves in A^-body simulations have been reported with re- 
markable coi ncidences for the radii of the main resonances. 

Puchs et al. (2005-) developed a similar argument for 
waves in the sheared sheet. They found that amplified 
trailing waves could excite fresh leading waves in their 
second-order theory. They proposed this mechanism as 
an alternative source of the amplitude excess noted above 
in the simulations by Toomrc & Kalnais (1991), and they 
speculated that the same mechanism may also account for 

* ICuzzi et al\ I I2010I) found evidence for similar behavior within Sat- 
urn's A ring. 



the larger than expected amplitudes of spirals in global A^- 
body simulations (see i )7.5[) . As both this mechanism, and 
that discussed in the previous paragraph, employ terms 
that are second-order in the perturbation amplitude, they 
become import a nt onl y when features are strong. 

Patsis et al. (|l99ll ) attempted to construct, by orbit 
superposition, self-consistent steady spiral waves of finite 
amplitude to match the observed non-axisymmetric pat- 
terns in specific galaxies. They experienced great diffi- 
culty in finding solutions near CR, and suggested that ei- 
ther this resonance or the 4:1 resonance^ marks the outer 
radius of the spiral. Their finding is not unexpected for 
two reasons: (a) the dispersion relation for steady, tightly- 
wrapped, small-amplitude waves (Fig. |4]) predicts a for- 
bidden region around CR for all Q > 1, and (b) the non- 
linear dynamics of orbits in barred potentials (see ^9.3p 
finds only chaotic or bits near CR , whic h are unfavorable 
to self-consistency as lPatsis et al\ ()l99lt ) found. It should 
be noted that the difficulty of finding a self-consistent solu- 
tion near CR is a direct consequence of their assumption of 
a steady wave pattern; transient waves do not suffer from 
this problem (Fig. [5]) and rapidly growing groove modes 
(j^H) have peak ov e rdens ities a t CR. 

Tsoutsis et al. (|2009t ) and lAthanassoula et al\ (|2009( ) 
suggested that spirals in barred galaxies could be created 
by the slow migration of stars along Lyapunov manifold 
tube orbits emanating from the unstable Lagrange points 
at the ends of the bar{3 The ambitious hope of these 
preliminary papers is for an ultimate unified picture for 
the co-cxistence of bars, spirals and rings, all having the 
same pattern speed. 

7.5. Spirals in Global N-body Simulations 

Sellwood & Carlberg digsl and iSellwoodI (Il989bn re- 
ported that their global simulations manifested more vig- 
orous spiral activity than was consistent with amplified 
shot noise. A brief summary of some further results to 
support this claim is given here, and will be described 
more fully elsewhe re. 

As noted in ^14. 11 iToomrd ()1981h predicted the half-mass 
Mestel disk to be globally stable to small-amplitude dis- 
turbances. Thus A^-body simulations of this disk should 
exhibit no activity in excess of the inevitable swing-amp- 
lified shot noise. Fig. [TU] reveals that this is not the case. 
The ordinate shows the largest value oi 6 — AS/S from 
m = 2 disturbances in a sequence of simulations with in- 
creasing numbers of particles. The unit of time is Ri/Vc, 
where VcRi is the center of the inner angular momentum 
cut out. Thus the orbit period at this small radius is 27r. 

At t = 0, (5 cx N~^/^, as appropriate for shot noise, and 
swing-amplification causes an almost immediate jump by 
a factor of a few for all A^. When A^ = 5 x 10^, ampli- 
fied noise causes ^ 20% overdensities almost immediately. 
For larger A" the amplitude eventually rises to similar val- 

^ The resonance condition (eq. [SJl for a pure cos{m9) potential vari- 
ation also implies frequency commensurabilities at multiples of m. 
The small denominators that characterize the principal resonances 
(BT08, §3.3.3) arise at these "ultraharmonic resonances" only for 
non-circular orbits. A new family of orbits appears at the 4:1 reso- 
nance tha t closes after 4 radial oscillations (s ee S|9.3|and Sellwood & 
Wilkinson Tl993IV 

Interestingly, this idea har ks back to t he old "g arden sprink ler" 
model for spirals proposed bv lJeani II1923I) (see also lJean£|[l929l . fig. 
55 & pp 357-360). 
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Fig. 10. — The time evolution of the peak overdensity in a series 
of simulations of the half-mass Mestel disk with differ ent numbers of 
particles. The model was predicted bv lToomrd 1)19811 ) to be globally 
stable. The ordinate shows the maximum values of 5 on grid rings. 
The number of particles in each simulation rises by a factor 10 from 
AT = 5 X 10'' to iV = 5 X 10*, and the initial ampUtude is oc Af-^/^. 
Linear theory predicts the later amplitudes should have the same 
scaling. 



ues in later evolution, once the inner disk has developed 
a pronounced bar. But for the largest two values of N 
shown, a period of slow growth occurs after the initial 
swing-amplified surge, offering tentative support for the 
linear theory prediction of global stability, with the slow 
growth perhaps arising from the gra dual development of 
partic le correlations as described by iToomre fc Kalnaii 
(|1991h . Even in these cases, however, the amplitude rises 
more rapidly once 6 > 2%. 

Thus even in these highly restricted simulations, spiral 
activity always exhibits runaway growth - albeit more and 
more delayed as N is increased - behavior that is quite 
clearly not in accord with linear theory predictions. The 
rapid growth once 6 > 2% suggests that the behavior has 
already become non-linear in some respect at this modest 
amplitude. Note that the largest number of particles, N — 
5 X 10*, is within a factor of 100 of the number of stars in 
a real galaxy disk, where in reality the mass distribution 
is far less smooth, owing to the existence of star clusters 
and giant gas clouds. 

7.6. A Recurrent Instability Cycle? 

Sellwood & Lin (IIQSQD reported evidence for a recurrent 
instability cycle in their simulations of a very low-mass 
disk. They observed that each spiral pattern created sub- 
stantial changes to the distribution of particles at the Lind- 
blad resonances, which they suspected created conditions 
for a new global instability of the groove mode kind (see 
They explicitly demonstrated that the later features 
were indeed true instabilities, since when they continued 
a parallel simulation after randomly shuffling the particles 
in azimuth at some moment, the pattern speed of the next 
mode to appear was the same as that in the original sim- 
ulation that had not been interruptedEB Thus, each co- 

This behavior is inconsistent with the non-linear mode coupling 



herent wave leaves behind an altered DF that apparently 
provokes a new instability. 

The runaway growth in the larger N models shown in 
Fig. (TU] is not caused by a single unstable mode, but ap- 
pears to be a succession of separate coherent waves, each 
having a lower rotat ion rate and rea ching a higher ampli- 
tude than the last. ISellwoodI (|2000l ) reported similar be- 
havior in lower- A'^ simulations, and demonstrated that one 
of the waves did indeed cause strong scattering at the ILR. 
It should be noted that resonance scattering is a second- 
order effect, but the evidence shown in Fig. [TU] suggests 
that it becomes important at a relative overdensity of just 
AS/S ^ 2%. Exactly how the demise of one mode cre- 
ates the conditions for the next instability remains unclear, 
however. 

Since the only evidence for this behavior came from 
a (well-tested) A^-body code, it seemed best not to pur- 
sue th e idea fur t her u ntil supporting evidence could be 
found. ISellwoodI (|l994f ) therefore expressed the hope that 
evidence of resonance scattering could be found in the 
HIPPARCOS measurements of the local stellar kinemat- 
ics. The publication of the GCS with distances and full 
phase-space m otions of ^ 14 00 F & G dwarf stars (see 
Fig.[T]) enabled [SellwoodI (|2010D to show that the so-called 
Hyades stream is in fact caused by scattering at a recent 
ILR. 

This very recent evidence supports the idea that spi- 
rals in the local Milky Way, and presumably elsewhere, do 
in fact behave as the simulations had indicated. Further 
work is required to expose the details of the recurrence 
mechanism. However, it now seems misguided to search 
for a spiral instability as some devious sort of cavity mode 
in a smooth disk; i.e. the assumption of a featureless DF 
may have thrown the spiral baby out with the bathwater! 

8. BUCKLING INSTABILITIES AND WARPS 

The optically visible parts of galactic disks are usually 
remarkably thin and flat, whereas the more extended H I 
disks of many edge-on galax ies appear no ticeably warped 
with an integral sign shape (ISancis W97B- Stellar warps 



with an integral sign shape (laancisil 
(|Reshetnikov et aZj 120021 : ISaha et al. 



2009f ) are much less 



pronounced than the warps in the extended gaseous disks. 
The long-known w arp in the H I layer of the Milky Way 
(lOort et a/.lll958ll has been most recently analyzed by 
iLevine et al\ ( 2006[ ) , while the dust and stars of the outer 
disk appear to be distorted in a similar, though less exten- 
sive, shape (jRevle et a/.ll2009l ). Both gaseous and stellar 
warps are frequently asymmetric, as appears to be the case 
for the Milky Way. The fact that stellar warps usually fol- 
low the same warped surface as do the gaseous ones (see 
also I Cox et oil Il996f ) , is strong evidence that warps are 
principally a gravitational phenomenon. 

Warps are extremely common. H I observations of edge- 
on galaxies ([Bosmal Il99ll : iGarcia-Ruiz et al\ 1200211 ) find 
very high fractions of warps, and the true fraction must 
be even higher, since warps directed close to our line of 
sight may be missed. Warps can also be detected kinemat- 
ically even when the system is not edge-on. Their ubiquity 
suggests that warps are either repeatedly regenerated or 
long-lived. 

ideas discussed in ^7.41 
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Briggs (|1990( ) studied a sample of 12 warped galax- 
ies with high-quality 21-cin data, and found that galactic 
warps obey three fairly simple rules: 

1. The H I layer typically is coplanar inside radius R25, 
the radius where the B-band surface brightness is 
25 mag arcsec^'^ = 25/iB, and the warp develops 
between R25 and i?26.5 (the Holmberg radius). 

2. The line of nodes (LoN) is roughly straight inside 

^26.5- 

3. The LoN takes the form of a loosely-wound leading 
spiral outside i?26.5- 

Kahn & Woltjer (|l959f ) first drew attention to the wind- 
ing problem presented by warps (see also BT08 §6.6.1), 
and argued that while self-gravity would slow the differ- 
ential precession, it could not be stro ng enough to create 
a long-lived warp. iLvnden-Belll (|l965l ). on the other hand, 
suggested that warps could result from a persisting mis- 
alignment between the spin axis and the disk normal, i.e. 
a long-lived mode. 

This section describes the theory of bending waves in 
general before addressing warps more directly. Unsta- 
ble bending modes will be denoted buckling modes, al- 
though other names that come from plasma physics, such 
as hose, firehose, or hosepipe instabilities, are commonly 
used. 

8.1. Local Theory of Bending Waves 

Toomre (|1966D considered the bending stability of an 
infinite, thin slab of stars having a velocity dispersion ax 
in the x-direction and some characteristic thickness zq. 
The self-gravity of the slab causes it to bend coherently 
provided the vertical oscillations of stars are adiabatically 
invariant as they move over the bend, which requires the 
slab to be thin compared with the wavelength of the bend. 
Furthermore, if the bending amplitude is small, its effect 
on the horizontal motion is negligible. 

Toomre derived the dispersion relation for small-amp- 
litude, long-wave {kzQ <^ 1) distortions of the form 
h{x,t) ^ i7e*(fe^-"t): 

Lj^ ^2TTG^k\-(jlk'^, (13) 

where S is the vertically-integrated surface density of the 
slab. The first term on the RHS is the restoring force 
from the perturbed gravity caused by the bend, while the 
second is the inertia term due to the vertical accelera- 
tion needed for stars to follow the corrugations. The in- 
ertia term is destabilizing and outweighs the gravitational 
restoring force when A < Aj = a^/GT,, causing the dis- 
tortion to grow exponentially. The unstable range of the 
buckling instability is precisely complementary to that of 
the Jeans instability in the plane in the absence of rota- 
tion, which is unstable for wavelengths A > Aj (jToomrd 
I1966D . 

The dispersion relation (|13p predicts a buckling insta- 
bility, at sufficiently short wavelengths, for any razor-thin 
system. However, it does not hold for wavelengths shorter 
than or comparable to the actual vertical thickness of the 
slab. Since zq ~ (t^/(GS), where cr^ is the vertical veloc- 
ity dispersion, one expects that bending at all wavelengths 



will b e suppres s ed wh en Uzjox exceeds some critical value, 
which iToomrel (I1966D estimated to be 0.300 
'^^r^i (|l985l) carried through a linear normal mo de anal- 
ysis of the infinite, isothermal slab (Spitzcr 1942j; ICamml 
|1950D . which has the properties p = po sech^(z/2zo), = 
(j1/{2'kGYj), and S = Ap^zo. He assumed a Gaussian dis- 
tribution of velocities, with ^ cr^, and determined 
the range of unstable wavelengths as the slab was made 
thicker. He showed that the buckling instability could be 
avoided at all wavelengths when CTz > 0.293f7a;, in good 
agreement with Toomre's earlier estimate. At the marginal 
stability threshold, the last unstable mode has a wave- 
length ~ 1.2Aj. 

Galaxy disks are not, of course, infinite slabs subject to 
plane-wave distortions, but the radial velocity dispersion, 
(jfl,, which is larger than the azimuthal dispersion, could 
perhaps be great enough to drive a buckling instability. 
As the observed ratio of velocity disper sions for stars in 
the So lar neighborho od is (Tz/cfr ~ 0.6 ([Nordstrom et al\ 
20041). Pfoomrd (|1966( ) concluded that at least this region 
of our Galaxy is "apparently well clear of this stability 
boundary." It should be noted that the approximate solu- 
tion for the potential, which requires kR 3> 1, is stretched 
in this case, since Aj 7kpc when gr = 40 km s~^ and 
S = 5OM0 pc-2. 

When the disk is embedded in some external potential, 
arising from a halo or the distant bulge of the galaxy say, 
the dispersion relation for short-wavelength waves in a thin 
slab (|13p acquires an additional stabilizing term 

= ul,^ + 2^GE|/c| - alk\ (14) 
where vl^^ = \d'^'^eKtl dz%=^ (BT08 eq. 6.114). Taking 
this additional factor into account further reinfo rces local 
stability and global, axisymmetric simulations (ISellwoodI 
[1996a) confirmed that Toomre's conclusion holds every- 
where in an axisymmetric, but otherwise plausible, model 
of the Milky Way disk. 

As for WKB spiral waves f i|3.2|) . the dispersion relation 
(jl4p can be generalized to tightly-wrapped non-axisym- 
metric bending waves simply by replacing lo with lo — mJl, 
with the angular rate of precession of the bending wave 
that has m-fold rotational symmetry being VLp = uj/m. 
It should be borne in mind that since observed warps in 
galaxies are very far from being tightly wound, analyses 
that make this approximation yield results that are at best 
only indicative of the dynamical behavior. 

Vertical resonances between the bending wave and the 
stars arise where m{Q,p — fi) = ft ^2.6p . although the 
meaning of v here depends on the context. BT08 (§6.6.1) 
considered only razor thin disks, for which the internal 
oscillation frequency J/int is infinite and the resonances oc- 
cur where m(f2p — fi) = ±fext- In disks of finite thick- 
ness, the stars have a natural internal vertical frequency 
vf^^ ~ AirGpQ (exact in the mid-plane of an infinite slab), 
and for vertical resonances the appropriate value of v = 
(r^i^t -I- j/gxt)^/2. Henceforth, will be used to mean ei- 
ther this total frequency in a thickened layer or v^xt for a 
razor-thin sheet. 

Eq. (ITi|) is satisfied for both ±m(57p — 51), implying two 
possible pattern speeds known as fast and slow waves. 

1^ IKulsrud et ali JTOTll) and lFridman fc Polvachenkol iITQSI) consid- 
ered the bending instability in a constant density slab of stars with 
sharp edges. 
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Because the gravity term raises the Doppler-shifted fre- 
quency above v, waves in a cold disk {aji — 0) extend 
away from the vertical resonances, and do not exist in the 
region between them that includes CR. The slow wave, 
which is a retrograde pattern for m = 1, is of most inter- 
est because the unforced precession rate, Hp = Vl — has 
much the milder radial variation. 

Much like density waves, neutrally stable bending waves 
propagate in a disk in t he radial direction at the group 
velocity (|Toomrelll983l ) 

du _ sgn(/c)7rGE — ka\ 

Since the denominator is negative for slow bending waves, 
trailing waves {k > 0) in a cold disk {aft, = 0) propagate 
inwards, while leading waves (fc < 0) propagate outwards. 
As for waves in non-uniform whips, the bending ampli- 
tude rises when the wave propagates into a region of lower 
surface density, and vice versa. 

Unfortunately, a full description of wave propagation in 
a sheet of finite thickness requires a solut ion of the lin- 
earized Boltzmann and Poisson equations (|Toomrelll966l : 
lA^aki 1985; Weinberg 1991; Toomrc 1995). The results are 
not analytic and surprisingly more complex than equation 
(|14p for the razor-thin case. Because the vertical potential 
of the disk is anharmonic, stars whose vertical oscillation 
takes them far from the mid-plane have lower vertical fre- 
quencies. Thus, any Doppler-shifted frequency m{flp — fl) 
of the bending wave will be in resonance with some stars 
that will damp the wave by converting wave energy into 
increased random motion. In general, short-wavelength 
mo des kzn > 0.5 damp i n less than one wavelength ( Wein- 
berg Il99ll : jTbOTnre 19951)" 

while long-wavelength modes 
kzo 'C 1 can propa g ate la rge distances. 

Sellwood et al. (|l998l ) studied the propagation and 
damping of axisymmetric bending waves in disks having 
both finite thickness and non-zero random velocities. 
Waves launched from the center of the disk at a fixed fre- 
quency propagated outwards and were damped as they 
approached the vertical resonance. As a result the disk 
thickened over a small radial range, with the peak occur- 
ring just interior to the resonance. 

8.2. Global Bending Modes 

Because of the complication caused by finite thickness, 
analytic work almost always adopts the razor-thin approx- 
imation, which also requires an = 0, since thin disks with 
velocity dispersion are buckling unstable^ With these 
assumptions, the disk can be approximated by a finite 
number of gravitationally coupled circular rings of matter, 
each having the appropriate angular momentum, which al- 
ways admits a discrete spectrum of bending modes. The 
real modes of the continuum disk can be distinguished by 
showing they are independent of the number of rings em- 
ployed. 

The gravitational restoring forces arc correctly captured 
in this approach, but the lack of random motion omits 
the additional coupling between adjacent mass elements 
caused by the epicyclic motions of the stars. This extra 

A razor thin disk is not destabilized by orbital motions with no 
velocity spread. 



mec hanism further stiffens the di sk's resistance to bend- 
ing (jPebattista fc SellwoodI [19991 ) . especially in the high 
density i nner regions where radial epicycles are larger. 

Hunter fc Toomre (|1969( ) developed the coupled ring 
approach to study the bending dynamics of rotationally- 
supported, razor-thin disks with no random motion and 
no halo. They were able to prove that general disks of 
this kind have no axisymmetric (m — 0) or warping (to = 
1) instabilities, but they could not extend their proofs to 
higher sectoral harmonics. In fact they noted that were the 
disk composed of two equal counter-rotating populations 
of stars (on circular orbits) it would be buckling unstable 
for all TO > 2. 

Hunter fc Toomre (119691 ) also studied the particular case 
of the sharp-edged Maclaurin disk in which all stars orbit 
at the same angular rate. When the sense of rotation is 
the same for all stars, the disk is stable to all bending 
waves, an d there is a simple s et of discrete neutral modes 
for all TO. iPolvachenkol (|l977f ) extended their analysis to 
Maclaurin disks with random motion (the Kalnajs disks) 
and was able to solve for the complete spectrum of bending 
modes. Needless to say, the introduction of random veloc- 
ities in this razor-thin system caused buckling instabilities 
to appear for all to. 

Since t he Maclaurin disk b e ars lit tle resemblance to real 
galaxies, 'Hu nter fc Toomri (|1969') modified the disk to 
blur the sharp edge. They demonstrated that discrete 
warp modes in a cold, razor-thin disk can exist only when 
the edge is unrealistically sharp. Note that all isolated 
disks admit two trivial zero-frequency modes: a vertical 
displacement of the entire disk and a tilt of the disk 
plane to its original direction. More interesting standing 
wave solutions (modes) require traveling waves to reflect 
off the disk edge, but a realistic disk with a fuzzy edge does 
not reflect bending waves, because the group velocity de- 
creases with the disk surface density (eq. [T5|) an d a wave 
packe t in a cold disk will never reach the edge ([Toomrd 
[l983l ). 

8.3. Simulations of Buckling Modes 

Sellwood & Merritt (|1994[ ) used iV-body simulations to 
study the global instabilities of hot disks with no net ro- 
tation, i.e. with half the particles counter- rotating. (See 
for buckling modes of more normal disks with large 
net rotation.) The form and vigor of the principal insta- 
bilities in any one of their models varied with the bal- 
ance between radial and azimuthal pressure and with disk 
thickness: an in-plane lop-sided instability was the most 
disruptive for cool models (SJH]). The radially hotter thin 
disks were disrupted by axisymmetric bending instabili- 
ties (bell modes), as illustrated in Fig.[TTJ The instability 
creates a thick inne r disk resembling a pseudo-bulge ( Kor- 
mendy & Kennicutt l2004D . Verv thin disks also buckled in 
an TO = 2 "saddle" mode, and an to = 1 warp instability 
was detectable in some models, but never dominant. 

Remarkably, instabilities in a counter-streaming model 
having intermediate radial pressure caused rather mild 
changes and led to an apparently stable, moderately thin, 
and almost axisymmetric disk. The in -plane velocities in 
this model resemble those reported by (jRubin et aLlll993 ) 
for the SO galaxy NGC 4550, and indicate this galaxy could 
be stable even without large quantities of dark matter. 
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Fig. 11. — The development of a buckling instability in a simple model of an isolated stellar disk with Q ~ 2. The disk is KT/4 model 
described by Sellwood & Merritt (1994), with equal numbers of particles orbiting in each direction, but here the simulation uses a code of 
much higher spatial resolution. The buckling mode is axisymmetric and, while non-axisymmetric features were permitted by the code, none 
developed. The disk mass and radial scale are unity, and G = 1; the orbit period at i? = 2 is 16 in these units. 



The stability of this end product also demonstrated that 
thin axisymmetric systems with modest radial pressure are 
more stable to bending modes than those having isotropic 
or radially biased DFs . 

Merritt & Sellwood ()1994[ ) proposed a criterion for the 
stability of a stellar system to buckling modes that can 
be applied globally. A particle moving at speed i; in a 
mildly bent sheet with characteristic wavenumber k expe- 
riences a periodic vertical forcing at frequency ki and, like 
any harmonic oscillator, the phase of its response depends 
upon whether the forcing frequency is greater or less than 
its natural vertical frequency. If kx < v, the density re- 
sponse of the system to an imposed perturbation will be 
supportive, and the disturbance can be sustained or even 
grow. However, ii kx > v for most particles, the overall 
density response to the perturbation will produce a poten- 
tial opposite to that imposed and the disturbance will be 
damped. Thus short waves are stable. They showed that 
their proposal successfully accounted for their TV-body re- 
sults, for the behavior of the infinite slab at short wave- 
lengths, and for the apparent absence of elliptical galaxies 
more flattened than E7. 

Sellwood ()1996bf ) found a long-lived bending oscillation 
in an A'^-body simulation of a warm disk that was con- 
strained to be axisymmetric. Apparently the system was 
able to support a standing wave, at most mildly damped, 
between the center and the edge of the disk, even though 
the stars had random motion in a disk of finite thickness. 
The frequency of the bending or flapping mode was low 
enough to avoid vertical resonances with almost all disk 
particles , in a greement with the requirement stated by 
iMathuij (|l99Cll ). While this result remains an isolated cu- 
riosity, since the flapping mode would probably be quickly 
damped in a halo, it is a clea r counter-example to the ar- 
gument bv lHunter k, Tooinri (|1969. ) that realistic disks do 
not possess global bending modes, which apparently holds 
only for disks without random motion. 

8.4. Disks in Halos 



Since the paper by iHunter fc Toomr3 (jl969t ) , ideas of 
warp formation have relied in some way on the interaction 
bet ween the disk and its dark matter halo. Dekel & Shlos- 
man (|1983D and lToomrd (119831 7 suggested that a flattened 
halo misaligned with the disk co uld form a long-las ting 
warp. ISparke fc Casertanol ()1988f ) andlKuijkeij (|1991h ob- 
tained long-lived warps (dubbed modified tilt modes) 
of disks in rigid, misaligned h alos, which wer e insensitive 
to the details of the disk edge. iLovelacj (|1998[ ) studied the 
tilting dynamics of a set of rings also in a rigid halo, but 
assumed that the inner disk lay in the symmetry plane of 
the spheroidal halo. 

Dubinski & Chakrabarty (120091 ) noted that the dark 
matter halos that result from cosmic structure formation 
simulations are usually aspherical, with frequent misalign- 
ments between the principal axes of the inner and outer 
halo. The disks in their simulations warped nicely when 
forced with slowly rotating, but otherwise rigid, perturb- 
ing fields representative of such halos. 

However, dark matter halos are not rigid, and a re- 
sponsive halo al ters the dvnamics in several wavs. Nel- 
son & Tremaine"ll995D showed that were the inner disk 
misaligned with the principal plane of the fla ttened halo, 
as supposed by S parke fc Casertanol ([1988) , its preces- 
sion would be damped through dynamical friction, bring- 
ing the disk into alignment with the halo on time scales 
much shorter than a Hubble time. But a more compelling 
objection to the modified tilt mode emerged from i V -body 
simulations with live halos: iDubinski fc KuiikenI (Il995f) 
found that the warp did not survive while iBinnev et al\ 
(199^ showed that the inner halo quickly aligns itself with 
the disk, not vice versa. The large store of angular momen- 
tum in the disk maintains its spin axis, but the pressure 
supported inner halo can readily adjust its shape slightly 
to align itself with the disk. 

8.5. Misaligned Infall 

The idea that galaxy warps are manifestations of eternal 
warp modes seems doomed by the damping effect of a live 
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halo. But slowly evolving warps remain viable, provided 
that suitable external perturbations occur in enough cases. 

In hierarchical galaxy formation scenarios, late infalling 
material probably has an ang ular momentum ax i s mis- 
aligned with the disk spin axis. lOstriker fc BinnevI (|1989f ) 
therefore proposed that warps arise due to the slewing of 
the galactic potential as material with misaligned angular 
m omentum is accreted. S tructure formation simulations 
bv lQuinn fc BinnevI ()1992l ) confirmed that the mean spin 
axis of a galaxy must slew as late arriving material rains 
down on the early disk. The less-than-critical matter den- 
sity in modern ACDM universe models implies that infall 
is less pervasive at later times, but it manifestly contin- 
ues to the present day in gravitationally bound environ- 
ments llSancisi et aLll2008l . and chapter by van Woerden & 
Bakker). 

Jiang feBinney dTool and lShen fc SellwoodI (l2006h pre- 
sented results of experiments in which a disk was sub- 
jected to the torque from a misaligned, massive torus at a 
large radius. This well-defined perturbation is a very crude 
model of an outer halo that is rotationally fiattened, and 
having its spin axis misaligned with that of the disk. It 
is misaligned and farther out because, in hierarchical sce- 
narios, the mean angular momentum of the later arriving 
outer halo is probably greater and misaligned from that of 
the original inner halo and disk. The accretion axis is, in 
reality, likely to slew continuously over time, so a model 
with a constant inclination is som ewhat unrealistic . 

Rather than striving for reahsm. lShen fc SellwoodI ()2006( ) 
used this simple forcing to reach an understanding of how 
the warp develops and why the LoN usually forms a loosely- 
wound, leading spiral. The inner disk maintained a coher- 
ent plane because of both self-gravity and random mo- 
tion. The torque arising from the misaligned outer torus 
caused the inner disk, and the aligned inner halo, to pre- 
cess rigidly even though the torque increased with radius, 
but the outer disk beyond ~ ARd started to warp. 

As soon as the outer disk became misaligned with the 
inner disk, the strongest torque on the outer parts of the 
disk arose from the inner disk. The torque from the inte- 
rior mass was responsible for the leading spiral of the line 
of nodes, even though the adopted external field would 
have produced a trailing spiral. The fact that the LoN of 
most warps forms a leading spiral over an extended radial 
range seems to imply massive disks. 

Even though the disk precessed due to the external 
torque, its motion was barely damp ed over many Gyr, 
in con trast to the expectations from iNelson fc Tremaind 
(|1995[ ). Damping was weak because the slow precession 
rate allowed the inner halo to remain closely aligned with 
the disk, which therefore caused little drag. The weak 
damping seemed to be caused more by the relative pre- 
cession of the inner and outer parts of the halo. Also the 
warp evolved slowly as the layer settled to the main plane 
at gradually increasing radii, in apparent agreement with 
the decreasing outward group velocity. 

Fi g. [T^ shows the H I observation of NGC 4013 bv Bot- 
tema (|1996[). together with the wa rp obtained in the sim- 
ulation bv IShen fc SellwoodI (|2006[ ). Their simulations re- 
vealed that the warp persisted for cosmologically inter- 
esting times, even when the external forcing field was re- 
moved. Thus the persistence of warps is not nearly as per- 




FlG. 12. — The up per panel shows the observed H I warp of NGC 
4013 (BottcmaT99D), and the lower panel the warp in the simulation 
by Shcn & ScUwood (2006) that closely resembles it. The length unit 
on the axes is the scale length of the exponential disc. 



plexing as previous studies had suggested. Furthermore, 
the model had a flat inner disk and the warp in the outer 
disk matched all of Briggs's rules quite well. 

A fixed outer torus is clearly unrealistic and the halo axis 
probably shifts continuo usly or episodically, as argued by 

uinn fc BinnevI (|l992l ). making warp lifetimes a side is- 
sue. Warps formed this way can be repeatedly regenerated 
when a new infall event happens. Since cosmic infall and 
mergers are more likely to happen in a denser environment, 
warps can be induced more frequently in such an environ- 
ment . which is consistent with the statistics ( Garcia- Ruiz 
et aL .2002bl ). 

8.6. Warps Driven by Tides 

Tidal interactions between galaxies could be an addi- 
tional cause of warps in disks, and one that is quite likely 
to produce asymmetrical warps. This idea has been ex- 
plored most fully to explain the warp of the Milky Way's 
(MW) disk that results from the proximity of the Magel- 
lanic Clouds, especially the Large Cloud (L MC). The ori- 
entation of the principal axes of the warp (iLevine et al. 
2006 ,) at least seems favorable to this hvpothesis. iBailin 
l|2003f ) suggested that the Sagittarius dwarf galaxy is an- 



other possible culpri t. 

Hunter fc Toomre (119691 ) , in models that did not include 
dark matter, concluded that the LMC could be responsible 
for the warp of the MW, but only if it was a lot more mas- 
sive than was then suspe cted and had recen t ly pass ed close 
to the edge of the disk. iGarcia-Ruiz et al. I (I2002ar i tested 
the hypothesis in fully self-consistent simulations that in- 
cluded live halos and used updated information about the 
distance and motion of the LMC. They concluded that nei- 
ther the amplitude, nor the orientation, of the warp in the 
disk of the MW was consistent with the tidal hypothesis. 
iWeinberg fc Blitz, (,2006i) . on the other hand, suggested 
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that a co-operative response from the halo to the passage 
of the LMC could have generated the observed warp in the 
disk. Thus there is no consensus yet on the origin of the 
warp in the Milky Way. 



9. BARS 

Bars are common features in disk galaxies. An earlier 
review ( Sellwood fc Wilkinsonl[l993l ) of the vast topic of 
bars is now somewhat dated but, as it provides a still useful 
summary of the basic dynamics, the present article will 
update a few main points and the reader is referred to the 
earlier review for a more detailed discussion. 

9.1. Origin of Bars 

For a long time, dynamicists were struggling to under- 
stand the absence of bars in some disk galaxies f Ostriker 
& Peebles Il973f r but since reasonable models of luminous 
galaxies that include a dense bulge of moderate mass are 
now known to be stable O, the problem has become al- 
most the opposite! Strong bars are seen in many galaxies 
whose mass distributions now appear unfavorable to the 
dynamical bar instability, as evidenced by a nuclear gas 
ring (see !j9.4[) . which can form only if the center is dense 
enough to have inhibited bar formation by Toomre's mech- 
anism. 

However, the fact that the center is dense today does 
not require that it was dense when the bar formed; sec- 
ular inflow of gas ( ij9.4l) can build up the central density 
after the bar has formed. Alternatively, the bar c ould have 
grown in size and slowed tr hough disk evolu tion (ISellwood 



1981 1: iBerentzen et al\\200T[ ) or halo friction ([Athanassoula 
20021 . see , or have been tr iggered by large density 



fluctuati ons in the disk (ISellwood ,1989aiV by tidal inter- 
actio ns (jNoguchil 1 19871 : IBerentzen et aLll2004l: ICurir et al\ 



120061 ) ■ halo substructure (IR omano-Di'az et al! '2008bl) or 
a non-axisymmetric halo (Dubinski & Chakrabarty 200^. 
Any of these considerations, or that described in the next 
paragraph, could plausibly reconcile the existence of bars 
in these galaxies with Toomre's stabi lizing mechan i sm. 

In an elegant piece of dynamics, iLvnden-Belll (|l979f ) 
demonstrated that the inner parts of galaxies are regions 
where eccentric orbits have a tendency to align themselves, 
which allows a bar to grow slowly through orbit trapping. 
The region where a cooperative response to a mild per- 
turbation occurs is where the overall radial density profile 
of the galaxy flattens into a more uniform core. In this 
region, the radial variation of — k/2 has a maximum at 
some non-zero radius and an infinitesimal bar pattern can 
have a pattern speed that allows two ILRs. Lynden-Bell's 
aligning mechanism, which operates best on the more ec- 
centric orbits, requires flp ~ fi^ — rifl/2 of orbits in the 
aligning region. As this pattern speed is much lower than 
that expected from the global bar instability f H4.1[) . the 
aligning mechanism offers an additional route to bar for- 
mation in otherwise globally stable disks. Although the 
coope rative region has a small radial extent, iLvnden-Belll 
(Il979t) suspecte d the bar could be much larger. 

Erwin ('2005') presented a useful study of bar properties 
and pointed out that bars in late-type galaxies are often 
much smaller relative to the disk size than arc those formed 
in sumulations. 

9.2. Frequency of Bars 

Strong bar s are clearly visible in 25% — 30% of disk 
galaxies (e.g. iMasters et all 120100 . and the fraction rises 
to > 50% when more objective criteria arc applie d to red 
or near-IR images ([Eskridge et aL.,2000, : .Marinova fc Joged 
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20071: Reese et al]|2007l . and further references cited below) . 



Barazza et al.\ (|2008|) found a higher bar fractio n in later 
Hubble types, while iMendez-Abreu et all (|2010 ^ found no 
bars in either very luminous, or very faint, galaxies. 

Whatever may be the mechanism responsible for the for- 
mation of bars in real galaxies, none of the above sugges- 
tions makes a clear prediction f o r the frequency of bars. 
iBosmal (|l996l ). ICourteau et al\ |2003), and others have 
pointed out that barred galaxies seem little different from 
their unbarred cousins in most respects, e.g. they lie on the 
same TuUy-F isher relation. Minor syst ematic differences 
do exist: e.g. iDavoust fc Continil (|2004[ ) note that barred 
galaxies seem to have smaller mass fractions of neutral H I 
gas, but this seems more likely to be the result, rather 
than the cause, of the bar. The a nti-correlation of bar 
frequency with the bulge half-light (iBarazza et aLl[2008f ) 
possibly results from Toomre's stabilizing mechanism, but 
this cannot be the whole story because some near-bulgeless 
disks are unbarred while other barred disks have massive 
bulges. 

If no dynamical property, other than their eponymous 
one, can be identified that cleanly separates barred from 
unbarred galaxies, then the existence of a bar in a galaxy 
may possibly be d etermined by external f actors, such as a 
chance encounter. lElmegreen et al\ (|1990|) reported an in- 
creased fraction of bars in groups a nd cluster s, but more re- 
cent work (,Barazza et aLll20d9l : lLr at 20091 : [Aguerri et al\ 
I2009D has found little or no variation of bar fraction with 
environment. It is also possible that the bar fraction could 
be changi ng with time: different groups disagree ( Jogee et 
a/. T2004l : ISheth et a/.ll2008l ) . probablv because observations 
of galaxies at significant look-back times are subject to 
systematic difficulties due to band-shifting and c hanging 
spatial resolution (see also lElmegreen et a/.ll2007f ). 

A radical alternative is to regard bars as transient fea- 
tures that form and decay, and that the current fraction 
of bar red galaxies represents the duty cvcle ( Bournaud et 
al. 200^ Bars in A^-body simulations are dynamically 
rugged objects that appear to last indefinitely. Of course, 
they could be destroyed by a merger event, for example, 
although not much in the way of a cool disk would sur- 
vive such an ev ent. Norman and his co-workers ( Pfen- 
niger & Norman ll990l : iNorman et aLlll996l ) have argued 
that bars can be destroyed by the accumulation of mass 
at their centers, which may lead to a pseudo-bulge an d /or 
a hot inner disk. Howeve r, IShen fc SellwoodI (|2004D and 
lAthanassoula e^"al] (j2005 l) found that unreasonably large 
and/or dense mass concentrations were r equired to cause 
their b ars to dissolve. The simulations bv lBournaud et ol] 
(|2005f ) uniquely show that gas accretion may aid the disso- 
lution of the bar and, with star formation, could recreate 
a cool disk that would be needed to make a new bar. Even 
if this behavior can be confirmed by others, the model re- 
quires very substantial gas infall. Moreover, the continued 
existence of the hot old disk and the build-up of a dense 
center makes every cycle of this speculative picture harder 
to achieve. 

In the distant future, galaxy formation simulations may 
have the quality and resolution perhaps to be able to 
predict the correct bar fraction, and thereby reveal their 
cause. 



9.3. Structure of Bars 

This section gives a brief description of a few important 
aspects of the orbital beha vior in large-amplitud e bars , 
and the reader is referred to ISellwood fc WilkinsonI (|1993D 
for a more comprehensive discussion. Weak bars can be 
treated using epicycle theory (BT08 §3.3.3). Most early 
orbit studies in strongly barred potentials considered mo- 
tion confined to the plane perpendicular to the rotation 
axis. Even though 3D motion is much richer, the funda- 
mental structure of bars is most easily understood from 
in-plane orbits. 

Since bars are believed to be steadily-rotating, long-lived 
objects, it makes sense to discuss their structure in a frame 
that co-rotates with the potential well at the angular rate 
Op. A rotating frame has the effective potential 



\^IR\ (16) 
where $ is the potential in an inertial frame. 

The effective potential surface in the disk plane (eq. [T5|) 
resembles a volcano, with a central crater, a rim, and a 
steeply declining flank. The crater is elongated in the di- 
rection of the bar, and the rim has four Lagrange points: 
two maxima, Z/4 and L5, on the bar minor axis and two 
saddle points, Li and L2 on the bar major axis. (The fifth 
Lagrange point, L3, is the local potential minimum at the 
bar center.) Because of Poisson's equation, the density 
contours of the bar must be more elongated than those of 
the inner ^cff- 

Neither E nor is conserved in non-axisymmetric po- 
tentials, but Jacobi's invariant Ij (eq.lH]) is conserved even 
for strong bars that rotate steadily. Since Ij — i|t>p-|-$off, 
where v is the velocity in the rotating frame, contours of 
<I>cff bound the possible trajectories of stars having Ij less 
than that contour value. Stars that are confined to the 
bar, which also have Ij < <I'ofr(ii), are of most interested 
here. 

A periodic orbit is a possible path of a star in the 
rotating frame that retraces itself, usually after a single 
turn around the center, but always after a finite number of 
turns. Because the orbits close, the orbital period in the 
rotating frame is commensurable with the radial period 
and these orbits are also described as resonant orbits of 
the (strongly non-axisymmetric) potential. 

Periodic orbits can be either stable, in which case a star 
nearby in phase space oscillates (librates) about its parent 
periodic orbit in an epicyclic fashion, or they are unsta- 
ble, in which case the trajectory of a star nearby in phase 
space diverges exponentially (at first) from the periodic 
orbit. The orbits of stars that librate around a periodic 
orbit are known as regular orbits, those that do not li- 
brate about any periodic orbit are known as irregular or 
chaotic orbits. Chaotic orbits have only a single integral, 
Ij, while regular orbits have an additional integral (two 
more in 3D) that confines their motion to a hypersurface 
of smaller dimension in phase space. Regular orbits are 
the more interesting because the star's orbit can be more 
elongated than the potential surface that confines it, which 
is of great value when building a self-consistent bar model. 

The main features can be illustrated in the simple po- 
tential (c/. BT08 eq. 3.103) 
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Fig. 13. — Examples, in a rotating bar potential, of important 
periodic orbits that close after two radial oscillations for every turn 
about the center, the 2:1 resonant families. Those orbits elongated 
parallel to the bar axis (horizontal) are members of the xi family. 
The X2 orbits are elongated perpendicular to the bar. 



where = x'^ +y'^, Rc is a core radius inside of which the 
potential is approximately harmonic, (? < 1 is the flatten- 
ing, and fo is the circular speed at large R when q — 1. As 
in BT08, the values are: vo = I, q = 0.8, Rc — 0.03, and 
rip = 1. With these parameters, the major-axis Lagrange 
points lie at a distance Rl — 0.9996 from the bar center. 

The 2:1 resonant periodic orbits shown in Fig. [T^] have a 
range of Ij values, but all close in the rotating frame after 
two radial oscillations for every turn about the center. Or- 
bits of the main family, denoted xi, are elongated parallel 
to the bar (horizontal in the Figure), and are referred to 
as the backbone of the bar, since the majority are sta- 
ble. Stars on these orbits all move in the same direction as 
the bar rotates, but have shorter periods because they are 
interior to the Lagrange points (or more loosely, they are 
inside CR). Some Xi orbits are simple closed figures, but 
others have loops near the outer ends where the star pro- 
gresses around the center of the galaxy more s lowly than 
does the rot a ting f rame. ISparke fc SellwoodI (|1987( ) and 
IVoglis et al\ (|2007l ) found that a large fraction of parti- 
cles in an A'^-body bar librate around them, and the same 
behavior is expected for stars in real bars. 

Another important family of 2:1 orbits appears near the 
center of the bar, but is elongated perpendicular to the 
bar axis, as also shown in Fig. 1131 This family, denoted 
X2, is almost always stable and is present in many real- 
istic bar potentials. Since Lindblad resonances ( ^12.4^ are 
defined for infinitesimal perturbations to axisymmetric po- 
tentials, the label ILR is very loose usage in strong bars. 
It is true that the X2 family appears in barred poten- 
tials only, but not always, when the axisymmetric mass 
distribution and pattern speed admit one or more ILRs. 
Furthermore, X2 orbits orient themselves perpendicular, 
while the xi family is parallel, to the bar axis, which is 
directly analogous to the abrupt phase shift that occurs 
in the response of a forced harmonic oscillator as the forc- 
ing frea uency crosses its natural freauencv (cf. Sanders & 
Huntleyllll^ At finite amplitude, the near circular orbit 
sequence acquires gaps at the 2:1 resonance bifurcations, 
which become broader as the bar amplitude rises. Thus 
the radial extent of the X2 family shrinks as the strength of 
the bar is increased, and it may disappear entirely. Even if 



one is careful to say that the appearance of the X2 family 
is the generalization of the ILR to strong bars, the radius 
of this resonance is still badly defined because the orbits 
that appear inside it can be quite eccentric. 

As I J approaches the value of $cff(ii), the time to com- 
plete a full turn in the rotating frame lengthens and ad- 
ditional orbit families appear. Orbits that close after any 
number of radial oscillations can be found in principle, but 
of these only the 4:1 resonant orbits are of dynamical sig- 
nificance to bar structure. As the period lengthens, the 
proliferation of orbit families causes a precipitous decrease 
in the stable r egions around e ach parent and chaotic be- 
havior ensues ([Chi rikov 197!^ . The onset of chaos near 
CR led lContopouloa (1980. ) to expect the density of a self- 
consistent bar to drop steeply near the major axis La- 
grange points, leading to the rule that the length of a bar 
is limited by CR . This rule predicts that the parameter 
(|Elmegreenlll996l ) 



n = RL/aB > 1, 



(18) 



where is the semi-major axis of the bar. In principle, 
self-consistent bar dynamics could allow bars with TZ ^ 
1, although empirical bar pattern speed estimates f ^9.5p 
mostly find that CR is in fact only slightly beyond the end 
of the bar. 

There are many more in-plane orbit families, but few 
are of dynamical importance to the structure of the bar. 
See lSellwood &: WilkinsonI ()1993f ) for a fuller account. 

The extension to 3D allows for many more resonances 
between the in-plane motion and the vertical oscillations. 
While there is a rich varietv of behavior f Pfenniger & 
Friedh 119911: ISkokos et aLl l2002f). the backbone xi fam- 
ily from 2D continues to be the most important, but now 
with a "tree" of orbits also librating vertically. The new 
periodic orbits that appear in 3D have a similar projected 
shapes as the in-plane xi family, but they also librate ver- 
tically a small number of tim es over t he same period as 
the motion in the plane. Pats is et al\ ([2002) highlighted 
the orbit families that they found to be of importance for 
the "boxy" appearance of edge-on bars (see also ^y^. 

9.4. Gas Flow 

When pressure and magnetic forces can be neglected, 
any mild dissipation will drive gas to move on stable pe- 
riodic orbits. An organized streaming gas fiow pattern 
is expected wherever the simplest periodic orbits over a 
range of energies can be nested and intersect neither with 
neighboring orbits, nor with themselves. Shocks, where 
pressure ceases to be negligible, must occur in flows either 
where periodic orbits self-intersect, or where gas flows on 
two separate orbits cross. Fig. [13] shows that were gas 
to flow in that adopted bar potential, shocks would be in- 
evitable because many orbits self-intersect (the loops) and, 
in particular, X2 orbits cross the xi family. Thus shocks 
are a general feature of cool (low pressure) gas flows in 

barso 

Full hydrodynamic simulat ions f e.g. [Roberts et a?]|l979l : 
lAth anassou Ialll992l : lFiI^[T999l ) are needed to determine the 
flow pattern. Shocks are off'set to the leading sides of the 
bar major-axis in models having reasonable parameters. 

Shocksjnav^e_avoide^^heii£ressure & 
Gerhard lToOTl) ! 
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iPrendergastj (|1962| ) seems to have been the first to asso- 
ciate the dust lanes in bars with the locations of shocks in 
the gas. 

Shocks convert some kinetic energy of bulk motion in the 
gas into heat, which is radiated efficiently. Furthermore, 
the offset location of the shock causes the gas to spend 
more than half its orbit on the leading sides of the bar, 
where it is attracted backwards towards the bar major 
axis, causing it to lose angular momentum^ Thus gas in 
the bar region must settle a little deeper into the potential 
well on each passage through a shock, i.e. the bar drives 
gas inwards, the angular momentum it loses being given 
to the bar. 

The inflow stalls where gas settles onto the X2 orbit fam- 
ily, which is found in bar models that have dense centers, 
leading to a build up of gas. This behavi or can be assoc i- 
ated with nuclear rings of dense gas (Regan et aLll2003 ). 
which are often the sites of intense st ar formation also 
(|Maoz et a/.ll200ll : [Benedict et aLll2002l ). If this were the 
whole story, the gas could not be driven any further in- 
wards, but there is both observational evidence, in the 
fo rm of spiral dust lanes and star formation (e.g. C arollo et 
al. Il998l ). and some theoretical work (|Wadall2o'oi[ ) to sug- 
gest that self-gravity causes inflow to continue. However, 
the existence of high gas density and rapid star formation 
in the nuclear ring indicates that only a small fraction of 
the gas continues inwards. 

Modeling the gas flow in a specific galaxy allows one 
to determine two properties of the bar that are hard to 
constrain otherwise. Estimating the gravitational poten- 
tial o f the galaxy f r om a photometric image plus a dark 
halo, IWeiner et al\ (|2001l ) fitted for both the disk mass- 
to-light (M/L) and ftp in the galaxy NGC 4123. These 
authors also described the procedure in detail. Results for 
a nuinber of other ga laxies were reported by Perez et~d\ 
(|2004D . IWein"CTl (|2004[ ). and lZanmar Sanchez et al\ (|2008f ). 
although these last a uthors were unable to obtain an en- 
tirely satisfactory fit. iLindblad et al\ ()1996f ) fixed the M/L 
and fitted only for flp. All these fits preferred rapidly ro- 
tating bars in heavy disks. Perez (2008.) confirmed that 
the best fit parameters of M/L and ftp were the same for 
both 2D Eulerian (Godunov) and 3D Lagrangian (SPH) 
hydrodynamic methods. 

Even though the quadrupole field of a bar decays quickly 
with radius, it can be strong enough to drive a spiral shock 
in the gas of the outer disk , as origina l ly dem onstrated by 
ISanders fc HuntlevI (|1976D . ISchwarzl (l98lh showed that 
when gas is modeled as inelastic particles, it is driven out- 
ward s to form an outer ring (see review bv lButa fc CombesI 
[Ml). However, it is unclear that spirals in the outer disks 
of real barred galaxies are the responses to the bar, and 
they may owe more to self-cxcited str uctures than to ba r 
driving (33 Sellwood & Sparkc 198S: iButa aLll2009t ). 
In addition, the outer spiral response to an imposed bar 
is not steady in modern simulations, with the shapes of 
the driven arms cycling through a broad range. For these 
reasons, the gas flow models fitted to individual galaxies 
should rely primarily on the fit within the bar and pay 
little attention to the outer disk. 

Kranz et al. (j2003l ) tried a similar approach, but fltted 

The opposite happens in shocks outside CR, where the gas gains 
angular momentum from the bar. 




Fig. 14. — Su mmary of direc t pattern speed measurements for 
bars coUected by ICorsin 3 1I2OOI) . The diagonal line shows TZ = 1 
and the dotted line TZ = 1.4. The crosses mark values for which 
error bars are unavailable. 



a spiral pattern instead of a bar, which may yield unre- 
liable results for two reasons: (i) The lifetimes of spiral 
arms are believed to be short (SJ7]) leading to broader res- 
onances and stronger gas responses than would arise in 
simulations that assume a slowly evolving pattern, and 
(ii) the observed spirals could be the superposition of sev- 
eral features with different angular rotation rates. As bars 
undoubtedly last for longer than do spirals and dominate 
the non-axisymmetric potential, fits to bars in galaxies are 
likely to yield better disk mass estimates. 

9.5. Bar Pattern Speeds 

Tremaine & Weinberg (1984a') devised a method to mea- 
sure the pattern speed of a bar directly from observations 
of a tracer component, which must obey the equation of 
continuity. Their original method assumes that the galaxy 
has but a single pattern, and would yield a misleading re- 
sult were there more than one pattern, each rotating at a 
different angular rate. 

The stellar component of early-type barred galaxies is 
believed to obey the equation of continuity because these 
galaxies have little dust obscuration and no star forma- 
tion. They also rarely possess prominent spirals in the 
outer disk. Results of many studies using this method 
for ear ly-type barred galaxies were summarized bv lCorsinil 
(^2008") and are shown in Fig. [TH While some individual 
measurements a re quite uncertain, the da. t a seem to favor 
1 < 7^ < 1.4. IChemenin fc Hernandez (|2009t ) found a 
counter-example in a low- luminosity galaxy. 

Fathi et al. (120091) andlMeidt et a l. (2009) applied the 
method of lTremaine fc WeinbereT i 1984al) to ionized and 
molecular gas, respectively. Both groups argue that this 
is valid, even though the separate gas components do not 
obey the continui ty equation that undelies the method. 
iFathi al] (|200^ generally find fast bars. Meidt et al\ 
()2008f ) generalized the method to atte mpt to measure ra - 
dial variations in the pattern speed and lMeidt et al\ (|2009t ) 
found suggestions of pattern speeds that are lower at large 
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radii than those near the center. 

Other methods can yield indirect estimates of bar pat- 
tern speeds. Fits of models of the gas flow ( ii9.4l) have been 
report ed for a few galaxies, finding TZ ^ 1.2. Athanassoula 
(|1992D argued that the shapes and locations of dust lanes 
in bars also seem to suggest that 7?. ~ 1.2. Associating a 
ring in a barred gala xy with the location of a major res- 
onance with the bar (jButa fc CombejllQQGT ) yields, with 
kinematic informat i on, an estimate of the pattern speed. 

Rautiainen et al. ()2008l ) computed models of the stellar 
and gas (sticky particles) responses to forcing by photo- 
metric models of 38 barred galaxies, in which they assumed 
that the entire non-axisymmetric structure rotated at the 
same pattern speed. They attempted to match the model 
to the visual morphology of the galaxy, and found a range 
of values for TZ. However, in most cases where TZ ^ 1, 
the fit is dominated by the outer spiral, which may have a 
lower angular speed than does the bar. 

9.6. Bars Within Bars 



Erwin & Sparke (120021 ) and others have found inner sec- 
ondary bars within the inner parts of > 25% of large- 
scale or primary bars. They reported that the secondary 
bar has a length some ~ 12% of that of the primary bar 
and the deprojected angles between the principal axes of 
the two bars appeared to be randomly distributed, sug- 
gesting that the two bars may tumble at dif f ering rates. 
This infe rence was supported bylCorsini et al\ (|2003D . who 
used the iTremaine fc Weinberd (|l984a ) method to show 
that the two b ars in NGC 2 950 c ould not have the same 
rotation rates; iMacieiewskil (|2006f ) used their data to ar- 
gue that the secondary bar has a large retrograde pattern 
speed. 

The theoretical challenge presented by these facts is 
substantial, and progress to unders t and t he dynamics has 
been slow. iMacieiewski fc Sparkd (|2000[ ) studied the or- 
bital structure in a potential containing two nonaxisym- 
metric components rotating at differing rates. However, 
it is almost certain that the secondary bar can n either ro- 
tate at a uniform rate ([Louis fc GerhardI 119881 ) nor can 
it maintain the same shape at all relative phases to the 
primary. 

Friedli fc Martinet ()1993f ) argued tha t gas was essential 
to forming secondary bars (s ee also e.g. iHeller et aL|[2001l : 
lEnglmaier fc Shlosman|[2004[ ) . However, some of the colli- 
sion less simulations re p orted by[Rautiainen fc Salo (1999) 
and iRautiainen et all (j2002f ) manifested dynamically de- 
coupled inner structures when the inner disk had high or- 
bital frequencies due to a dense bulge. The structure was 
more spiral like in some models, but others appeared to 
show inner bars that rotated more rapidly than the main 
bar. 

Debattista fc Shen ()2OO70 created long-lived, double- 
barred galaxy models in coUisionless A'^-body simulations 
having dense inner disks, which they described as pseudo- 
bulges. They followed up with a more detailed studv f Shen 
fc Debattista |2009( ) that also made some predictions for 
observational tests. The secondary bars in their models 
indeed rotated at non-uniform rates, while their shape var- 
ied systematically with phase relative to that of the pri- 
mary. These models prove that coUisionless dynamics can 
support this behavior, but it is unclear that their initial 



conditions mimicked those that have given rise to double 
barred galaxies in nature. 

The possible consequence of gas inflow in these galaxies 
has attracted a lot of attention. Shlosman e t al. (1989,) 
speculated that bars within bars might lead to gas inflow 
over a wide dynamic range of scales, from global to the 
parsec scale where accretion onto a black hole might cause 
AGN activity. W hile inflows may have been observed {e.g. 
iHaan e< aIll2009D . understanding of gas flow in these non- 
stea dv potentials remains rather preliminarv ( Macieiewski 
et aL 120021 : IHeller et aLll2007l 'l. 

9.7. Buckling of Bars 

Combes fc Sanders ()1981l ) first reported that the bars 
in their 3D simulations were thicker than the disk from 
which they had formed, and had acquired a pronounced 
"boxy" shape when viewed edge-on. Boxy isophotes in 
edge-on disk galaxies are now believed to be an indicator 
of a bar, as is support ed by kinematic evidence in the ga s 
(|Merrifield fc Kuiiken|[r999 : Bureau fc Athanassoulall2005[ ) . 

The reason the bar thickened was explained bv Raha et 
al. (jl99l[ ). who showed that bars are subject to the buck- 
ling instability f ^iS.ip . The bar buckles because the forma- 
tion of the bar created a structure supported by elongated 
orbits that stream along the bar in the near radial direc- 
tion. Even though the ingoing and outgoing stars stream 
on different sides of the bar, the effective averaged an has 
risen as a result w ithout changing tr^. The simulation by 
iRaha et al\ (|l99lD revealed that the buckling instability 
produced a large amplitude arch just before it saturated, 
after which the bar became thicker. The energy to increase 
vertical motion in the bar appeared to have been released 
by the further concentration of mass towards the bar cen- 
ter (see alsol l200ii ). It is delightful that the evolution of 
one instability, the bar-forming mode, should create a new 
structure, the bar, that is itself unstable. 

Bars still thicken in more recent simulations with grids 
having higher spatial resolution {e.g. Fig. O, but do not 
seem to exhibit the spectacular arch reported bv Raha 
et al. (|l99ll ) unless bi-svmmetrv is enforced. Low spa- 
tial resolution or significant gravity softening (which are 
equivalent) weakens the restoring force in eq. (|13p and ar- 
tificially increases A j , which is the characteristic length for 
instability. Stronger gravity causes the preferred buckling 
modes to have shorter wavelength allowing, say, an up- 
ward arch on one side of the center and downward arch on 
the other. Enforcing bi-symmetry prevents the bar from 
bending in this anti-symmetric manner, and forces it to 
buckle through the single arch mode. 

9.8. Dynamical Friction on Bars 

Friction between a rotating b ar and a massiv e halo was 
first reported many years ago (jSellwoodl Il980l ). but the 
implications for dark matter halos have fueled a renewed 
intense study of the topic. 

9.8.1. Theory 

Dynamical friction ()Chandrasekhaijri943f ) is the retard- 
ing force experienced by a massive perturber as it moves 
through a background sea of low-mass particles. It arises, 
even in a perfectly coUisionless system, from the vector 
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Fig. 15. — The dimensionless acceleration function V defined in 
eq. HIQI I for the case of a Gaussian distribution of velocities among 
the background particles. The function is negative because the ac- 
celeration is directed oppositely to the velocity. 



sum of the impulses the perturber receives from the par- 
ticles as they are deflected by its gravitational field (see 
Appendix). Equivalently, friction can be viewed as the 
gravitational attraction on the perturber of the density 
excess, or wake, that d evelops behind it as it moves, as 
was nicely illustrated bv lMuldeil (|l983f ). 

Chandrasekhar's formula (BT08 eq. 8.6) for the accel- 
eration of a perturber of mass M moving at speed vm 
through a uniform background, density p, of non-inter- 
acting particles having an isotropic velocity distribution 
with a ID rms velocity spread a, may be written as 



dv 



M 



(19) 



dt (7^ 
The Coulomb logarithm is defined in the Appendix, and 
the dimensionless function V is drawn in Fig. 1151 for a 
Gaussian distribution of velocities; other velocity distribu- 
tions would yield a different functional form. Physically, 
the retarding acceleration must vanish when the perturber 
is at rest and it must also tend to zero when the perturber 
moves so fast that the background particles receive only 
small impulses and the feeble wake lies far downstream 
from the perturber. Friction is strongest when the speed 
of the perturber is somewhat greater than the rms speeds 
of the background particles. 

The simplifying assumptions in its derivation invalidate 
application of eq. ([T9l) to the physically more interesting 
problem of friction in a non-uniform medium in which the 
background particles are confined by a potential well and 
interact with the perturber repeatedly. 

Repeated encounters between the perturber and the 
background particl es require the more soph is ticated treat- 
ment presented in iTremaine fc Weinberd (| 1984b). who 
adopted a rotating potential perturbation in a gravitation- 
ally b ound spherical halo of test pa rticles. They showed, 
as did lLvnden-Bell fc Kalnaia (|1972[ ) for spiral waves, that 
lasting changes to the orbits of the halo particles appear 
to second order in the perturbing potential, and can occur 
only at resonances. They derived a daunting formula for 
the torque on the halo caused by the perturbation that 
sums the contributions from infinitely many resonances. 
The contribution at each resonance is proportional to the 



gradient of the DF near the phase-space location of the 
resonance, in a manner that is directly analogous to Lan- 
dau damping. I Weinberg (|1985f ) computed the surprisingly 
large torque expected on a rotatin g bar, and his conclusion 
was confirmed in restricted t ests ()Little fc Carlber j|1991l : 
|Hernauist fc Weinbe_rdll 992'). 

Weinberg fc Katz (|2007) pointed out that friction is 
dominated by a few important resonances. They estimated 
the widths of these resonances for a perturbation having 
constant pattern speed and finite amplitude, and argued 
that immense simulations would be needed to populate 
the resonance with sufficient particles to capture the cor- 
rect net response. However, the loss of angular momen- 
tum from the perturber causes its pattern speed to change, 
and the resulting time-dependence of the forcing frequency 
is a much more important factor in broadening the reso- 
nances; thus friction can in fact be captured correctly in 
simulatio ns having moderate numbers of particles (^ Sell- 
wood ,2008ar r Note that the pattern speed of an orbiting 
satellite rises as it loses angular momentum, while that of 
a bar usually decreases. 

Despite the complicated language of resonant dynamics, 
the upshot is simply that the perturber induces a wake- 
like response in the halo, as was beautifully illustrated by 
IWeinberg fc Katj (|2007l . their Fig. 1). As for the infi- 
nite medium, friction can be thought of more simply as 
the torqu e between the p erturber and the induced halo 
response. ISellwoodl ()2006l , his Fig. 2b) shows the lag angle 
between the forcing bar and the halo response, which is 
about 45° when friction is a maximum and gradually de- 
creases to zero as the bar slows, until eventually friction 
ceases when the halo res ponse co-rotates with the bar. 

Li n fc Tremaine (Il983l) for an orbiting satellite, and Sell- 
wood (200&) for a rotating bar, demonstrated that the 
frictional drag on the perturbation scales with the mass 
of the perturber, M, the halo density, p, and the halo ve- 
locity dispersion, a exactly as in eq. (jl9p . Furthermore, 
the dimensionless function that describes the dependence 
on the angular speed of the perturber shares the general 
properties with V{x) that it is negative (for reasonable 
non-rotating halos), and must — > as a; — > oo, and that 
it should be oc X as a; — >■ 0. Including self-gravity in the 
halo response causes a further slight increase in friction, 
but does not otherwise change the behavior. 

9.8.2. Halo Density Constraint 

Fully self- consistent simulations of bar formation in a 
live halo by iDebattista fc SellwoodI (|1998l . 2000) showed 
that strong bars are indeed slowed rapidly. The fact that 
observed bars appear not to have been slowed f ij9.5p may 
imply an upper bound to th e density of the dark r natter 
halo in barred disk galaxies. IValenzuela fc KlvpinI (j2003( ) 
claimed a counter-example of a bar that does not experi- 
ence much friction in a "cosmologically-motivated" halo, 
even though their result disagreed with all others for strong 
bars (jO'Neill fc Dubinski .2003. : Athanassoula 200a) and 
with theory! 

Investig ation of their anomalous result bv Sellwood fc 
Debattista (12006 ) revealed that friction can be avoided 
temporarily if the gradient of the DF at the most important 
resonance(s) has been fiattened by earlie r evolution, which 
they described as a metastable state. iLin fc Tremainj 
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(|1983[ ) reported similar behavior as a result of driving the 
perturber at constant frequency for a protracted pe riod. 
They showed, as d id Sellwood fc Debattistal ()20Q6f ) and 
IViUa- Vargas et al\ (|2009l ). that the full frictional drag re- 
sumes after some d elay, the duration of which se ems to 
vary stochastically (jSellwood fc Debattistal [20091 ). De- 
layed friction can happen only in simulations of disks in 
isolated, smooth halos, since any reasonable amount of 
halo substructure, or a tidal encounter, disturbs the del- 
icate metastable state of the halo, causing friction to ap- 
pear with its full force. Thus simulations that do not find 
strong fricti on from moderately dense halos (e.g. Klvpin 
et al. 2009) hav e simply not been run fo r long enough. 

While .Debattista fc SellwoodI (|2000D argued for near 
maximal disks, and their requirement for a low halo den- 
sity is in agreement with the disk masses derived from 
fitting bar flow models ( ii9.4p . their constraint on the halo 
density may be specific to their adopted halo models. Thus 
additional careful studies of other halo models seem war- 
ranted. 

9.8.3. Halo Density Reduction by Bars 

While a full discussion of processess that may lower the 
dark matter density in the centers of halos is outside this 
review, a brief mention of the effect of bar friction is ap- 
priate here. 

Weinberg & Katz ()2002[ ) argued that the transfer of an- 
gular momentum from the bar to the halo could reduce the 
central density of the dark matter halo by a substantial fac- 
tor. However , the p ossible density reduction is quite mod- 
est (iHollev-B ockcl mann et a/.ll2005l : iMcMillan fc DehnenI 
1 



[20051 : ISellwood i2008al ) because the disk has only a finite 
amount angular momentum to give to the halo. Further- 
more, as the disk loses angular momentum, its mass distri- 
bution contracts, and the deepening potential well further 
compresses the hal o, which actual l y overwhelms the slight 
density reduction ([SellwoodI |20Q3[ : IColm et aLll2006l ). 



10. SECULAR EVOLUTION WITHIN DISKS 

Most of the behavior discussed so far, such as the non- 
linear evolution of instabilities, causes changes to the host 
galaxy on a dynamical time scale. The broad topic of 
secular evolution in galaxies describes changes that occur 
more grad ually, such as the secular formati on of ps cudo- 
bulges fsee lKormendvl[T993l : lKormendy fc Kennicuttll2004l . 
for exce llent reviews), the formation of rings (e.g. Buta fc 
Combes ll996D . or dynamical friction between components 
f TO.Sp . The discussion in this section concerns processes 
that scatter disk stars only. 

It has long been realized that old stars in the solar neigh- 
borhood have larger peculiar motions relative to the lo- 
cal standard of rest (hereafter LSR) than do young stars 
{e.g. IWielen' 1977; Nordstrom et aLll2004l: Aumer & Bin- 
ney ,2009) . Postulating that older stars were born with 
larger random speeds, say in a thicker disk, is unattrac- 
tive because it makes the present epoch of low velocity 
dispersion special. Some mechanism to scatter stars must 
therefore be invoked to create the larger random speeds of 
older stars. 

The trends presented in Fig. 1161 use the revised ages 
([Holmberg et ajj I2007D assigned to the GCS stars (see 
[Sodcrbloni.201Q . for a review). Figs. 2, 3 fc 4 of Aumer fc 
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Fig. 16. — The dispersion of stars in all three components, and 
the total dispersion of the GCS sample of solar neighborhood stars 
using revised a ges IHolmberg et a/.lj200'iTl . The radial, azimuthal and 
vertical components are crjj , av, pnd aw respectively, and (t^^^ = 
^Ij + + ■ The fitted straight lines do not include the first 
three orjlast two points. 



Binnev ( [2009D show the dispersions estimated as a func- 
tion of color on the main sequence; while blue stars are 
necessarily young, red stars are expected to have a range 
of ages. The dispersions of the s upposed oldest stars esti- 
mated bv lHolmberg et aj] (l2007^ are no greater than those 
estimated bv lAumer fc Binneyl ( 20091 ) for their reddest 
stars, suggesting that the "older" bins in Fig. IT6l include 
st arg hayin g a wi de range of ages, as argued by Reid et 
aL" l[2007^ 7 It should be noted that the reported disper- 
sions are simple second moments of the perhaps complex 
velocity distributions (see Fig. [T|) . 

A reliable determination of the variation of dispersion as 
a function of time could provide another useful constraint 
on the sc attering mechanism (e.g. BT08 ?8.4). Quillen 
fc Garnett lloooh and lSeabroke fc Gilmord ([20071 ) argued 
that the dispersion may saturate for stars above a certain 
age; with a much o lder surge to account fo r the highest 
velocities. However. lAumer fc BinnevI ([20091 ) found better 
fits to the data with continuous acceleration, and deduced 
(J ~ t^-^^, with t being the current age of the stellar gen- 
eration, in tolerable agreement, in fact, with 0.38 for the 
logarithmic slope of utot in Fig- [HI 
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Three principal scattering agents have been discussed: 
dense gas chimps in the disk, massive black holes in the 
halo, and recurrent short-lived spirals. Note that the first 
two are essentially collisional processes that accelerate the 
relaxation rate (see the Appendix), while the changes 
caused by spirals can increase random motion without 
leading to a more relaxed DF. As none in isolation fits 
the data, a combination of spirals and gas clumps seems 
to be favored. Minor mergers and the effects of halo sub- 
structure are discussed in t jlll 

An orbiting mass clump induces a collective spiral wake 
in the surrounding dis k that enhances its ma ss and size 
by substantial factors (jJulian fc Toomr^ll966D . a compli- 
cation that is ignored in many studies of cloud scattering. 
Since molecular gas is mostly concentrated in spiral a rms 
(jNieten et aLll2006l ; iGratier et a?.|[20Tol ; lEfremovll20Tol ). it 
is probably futile to draw a sharp distinction between spi- 
ral arms and the wakes of dense gas clumps, and a correct 
treatment would be to c alculate the effects of t he com- 
bined star and gas disk. iBinnev fc LacevI (|l988f ) took a 
step in this direction, but a full calculation may remain 
unreachable for some time if one tries to include a self- 
consistent treatment of the formation and dispersal of the 
gas clumps: molecular gas concentrations probably grow 
in the converging gas flow into a spiral arm, and are sub- 
sequently dispersed by star formation. 

Treating spirals and mass clumps in the disk as distinct 
scattering agents may be justified, therefore, if the wakes 
of cloud complexes can be lumped with spirals into a sin- 
gle scattering agent that is distinct from the clouds that 
caused them. At the very least, this simplifying assump- 
tion separates the problem into tractable pieces. 

10.1. Heating by Spirals 

Lynden-Bell & Kalnajs (Il972f ) showed that stars are 
scattered by a slowly changing potential perturbation only 
near resonances. More precisely, a spiral potential that 
grows and decays adiabatically, i.e. on a time-scale long 
compared with the orbital and epicyclic periods, will not 
cause a lasting change to a star's E and L^. Wave-particle 
interactions become important near the resonances, where 
stars experience secular changes through "surfing" on the 
potential variations at CR or through a periodic forcing 
close to their epicyclic frequency at the Lindblad reso- 
nances. Either case produces a lasting change to a star's 
orbit. 

The width of a resonance, i.e. the range of orbit fre- 
quencies of stars that are strongly affected, depends only 
on the amplitude of the potential when the perturbation 
is long-lived. But perturbations of shorter lifetimes have 
a broader range of frequencies and more stars experience 
lasting changes. 

The discussion in H2.4I and Fig. [3] indicates that stars 
that lose (gain) near the ILR (OLR) move onto more 
eccentric orbits, which is the root cause of heating by spi- 
rals. Exchanges at CR move stars to new orbits also, but 
with no change to the energy of non-circular motion, as 
discussed in iil0.2l 

Significant heating by spiral waves over a large part of a 
disk requires them to b e transient; a quasi-s teady pattern, 
of the type envisaged bv lBertin fc LinI ( 1996D say, will cause 
localized heating at an exposed resonance, while stars else- 



where will mov e through the patter n with out otherwise 
being a ffected. iBarbanis fc Woltie^ (I1967D. Carlberg & 
Scllwood~ifl985D and IBinnev fc La cevI (Il988l) calculated 
the hea ting caused by transient spirals. Jenkins fc Binnev 
(Il990j), iDe Simone et al\ (|2004r i. and iMinchev fc Quillen 
(|2006l ) presented numerical studies of the consequences for 



a disk of test particles subject to some assumed set of spi- 
ral wave perturbations. 

It is important to realize that the vertical oscillations of 
star s are little affec ted by spiral potential variations f H2.6l 
and lCarlberdll9 8?! ) . In the absence of heavy clumps that 
can redirect disk velocities f m0.3p . the increasing in-plane 
motions in simulations of initially cool, thin disks may ul- 
timately cause the velocity ellipsoid to become sufficiently 
anisotropic as to cause it to thicken through mild buck- 
ling instabilities fi jS.lD . This may account for claims {e.g. 
[McMillan fc Dchnen,2007j ) that disks thicken due to spiral 
heating. 

10.2. Churning by Transient Spirals 

Studies of th e metallicities and ages of nearby stars ( Ed- 
vardsson et aL Il993l : iNordstrom et al\\2004 ) found that 
older stars tend to have lower metalliciti es on average. As 
the ages of individual stars are disputed l|Reid et aLll2007l : 
iHolmberg et aLl [20071 ). the precise form of the relation is 
unclear. However, there seems to be general agreement 
that there is a spread of metallici ties at each age, which 
is also supp orted bv other studies (IChen et a/.ll2003l: Hav- 
wood 120081 : IStanghellini fc Havwoodll2010l ). The spread 
seems to be more than twice that expected from simple 
blurring of the gradient by stellar epicyclic excursions. In 
the absence of radial mixing, a metallicity spread amongst 
coeval stars is inconsistent with a simple chemical evolu- 
tion model in which the metallicity of the disk rises mono- 
tonically in each annular bin . 

Sellwood fc Binney (|2002D showed that scattering at CR 
causes very effective mixing. In a few Gyr, multiple tran- 
sient spirals caused stars to diffuse in radius. Churning 
of the stellar disk occurs at corotation of the spirals with 
no associated heating, and is able to account for the ap- 
parent metallicity spread with age. Roskar et al. (2008ab) 
presented more detailed simulations that included infall, 
star formation and fe edback that confirmed this behavior. 
ISchdnrich fc BinnevI ((2009) developed the first chemical 
evolution model for the MW disk to include radial churn- 
ing. 

10.3. Cloud Scattering 

Man y years before the discovery of gi ant molecular gas 
clouds, ISpitzer fc Schwarzschildl (|l953l ) postulated their 
existence to account for the secular heating of disk stars. 
iLacevI (|1984f ) extended their calculation to 3D and con- 
cluded that cloud scattering should cause the vertical dis- 
persion, Gz, to be intermediate between the radial, cr/j, and 
azimuthal, ti^, componentsl3 This result seems physi- 
cally plausible on energy equipartition grounds: scattering 
by massive clouds redirects the peculiar motions of stars 
through random angles, and therefore isotropizes the mo- 
tions as far as the epicycle gyrations allow. 

The ratio crfl/o"^ f» 2r2/K (BT08, eq. 8.117) is forced by epicyclic 
motions of disk stars. 
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Fig. 17. — The equilibrium axis ratio of the velocity elKpsoid of 
particles plotted as a function of the limiti ng range o f the p ertur- 
bation forces from the heavy particles. See ISellwoodI l l2008bl) for a 
description of the calculations. 



Despite the fact that redirecting pecuhar motions hap- 
pens much more rapidly than they can be increased by the 
same scatterers, the data do not reveal Lacey's expected 
axis ratio. The second moments of the velocity distribu- 
tion of solar neighbo rhood stars in the three orthogonal di- 
rections (Fig.[Ttland lWielen|[T977l :lA umer fc Binnevl[2009h 
satisfy the inequality cr^ < cr^ < a^. The ratio of the 
two in-plane components is in reasonable agreement with 
expectations from epicyclic motions, but the vertical com- 
ponent is the smallest, and this remains true for all groups 
when the stars are subdivided acc o rding to color or esti- 
mates of their ages. iGerssen et al\ (|200C ) also observed a 
flattened el li psoid in th e disk of NGC 2895 

Carlberg l|l987f ) and iJenkins fc BinnevI (|l990t ) there- 
fore developed the plausible argument that spirals drive 
up the in-plane components more rapidly than scattering 
is able to redirect those motions into the vertical direction, 
thereby account ing for the obser ved axis ratios of the ve- 
locity ellipsoid. ISellwoodI ()2000[ ) cited their argument as 
offering strong support for the transient spiral picture, but 
it now seems to be incorrect. 

Ida et al. (|l993l ) claimed that cloud scattering alone 
would lead to the vertical component being the smallest, 
with the precise axis ratio depending o n the local slope 
of the rotation curve. Their s imulations ( Shiidsuk e fc Ida 
and others ( e. a. IVillumsen,.1985. : Hanninen 



19991), and others (e.q. lVillumsen. 1985. : Hanninen fc Flvnn 
2002f). confi r med th eir expectation. 

Sellwood (j2008bD resolved this disagreement using sim- 



ulations of test particles in the sheared sheet (see 
Scattering by randomly distributed co-orbiting mass 
clu mps confirmed t he flattened velocity ellipsoid predicted 
bv llda et al\ (|1993[ ). Fig. [T7I reveals why Ida's prediction 
differs from Lacey's: Lacev. as others (Spitzer & Schwarz- 
schild7l953; Binncy & Laccy 1988), assumed that cloud 
scattering is local, but the In A term in the formulae in 
the Appendix implies that distant encounters dominate 
any scattering process in 3D, at least to distances a few 
times the disk thickness. Distant scatterers in the flat- 
tened geometry of a disk must predominantly affect the 
in-plane star velocities, and couple much less strongly to 



the vertical component. Fig. 1 171 shows the equilibrium ra- 
tio az/aji when scatterers beyond the finite range dmax 
artificially exert no forces. The ratio settles to something 
close to Lacey's energy equipartition prediction when none 
but the closest heavy clumps perturb the stars, but the 
equilibrium ellipsoid flattens in separate experiments as 
more distant clouds are included, tending towards Ida's 
result with no artificial cut off. 

Thus the shape of the local velocity ellipsoid. Fig. \W[ 
is apparently consistent with cloud scattering and does 
not, as seemed attractive, require concurrent spiral arm 
scattering. However, the data do not imply that spirals are 
unimportant: cloud scattering seems unable to generate 
the random speeds of th e oldest stars (e.g. iLacevI [l99TI : 
iHanninen fc Flvmil l2002f ) . and there are hints in Fig. [1^] 
of some evolution of the velocity ellipsoid shape that may 
demand a compound origin. 

10.4. Black Holes in the Halo 

The possibility that the dark matter halos of galaxies are 
made up of massive blac k holes (BHs) is not yet excluded. 
iLacev fc OstrikeJ (|l985f ) calculated the consequences for 
the stars in the disk of the MW, assuming the BHs to have 
orbits characteristic of a pressure-supported halo and to 
impart impulses to disk stars as they pass through the 
disk. The high speeds of their encounters with disk stars 
would cause the velocity dispersion to rise as ^1/2, while 
the predicted shape of the velocity ellipsoid is in reasonable 
agreement with that obser ved. 

Lacey & Ostriker ()1985f ) addressed a number of issues 
with their model, such as the X-ray accretion luminosity 
as the BHs pass through the gas disk, the accumulation of 
BHs in the galactic center through dynamical friction, and 
the survival of dwarf galaxies. They also acknowledged 
that it does not predict the correct variation of crtot with 
Galactic radius. BT08 (§7.4.4) added that the idea could 
be ruled out on the grounds that wide binary star systems 
would be disrupted too quickly. 

If disk scattering is dominated by spirals, as argued be- 
low, then scattering by BHs would be needed merely to 
redirect the peculiar motions into the vertical direction. 
This reviewer is not aware of any such calculation, but 
since passing BHs scatter stars in the direction perpen- 
dicular to their orbits, it seems unlikely that they could 
redirect peculiar motion without also increasing in-plane 
motions. However, if this expectation is too pessimistic 
and the desired axis ratio could be achieved with lower 
BH masses, many of the other objections are weakened. 

10.5. Discussion 

The extraor dinary p hase-sp ace structure of the solar 
neighborhood (jDehnenI 1998; N ordstrom et aZ.ll2004l . and 
Fig. [Ij indicates that there is little in the way of an un- 
derlying smooth component a nd the stellar dist ribution is 
broken into several "streams" (jBovv et a/.ll2009f ). The fea- 
tures are too substantial to have simply arisen from groups 
of sta rs that were born with similar kin ematics (e.g. Eggerij 
19961). as confirmed in detail ed studie s (iFamaev et aLll2007l : 
Bensbv et a0l2007l : lBovv"fc Hogg 2009^ 

Thus it is clear that the entire DF has been sculptured 
by dynamical processes. Were the large part of the spread 
in velocities caused by multiple scatterings off molecular 
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clouds in the disk, or off black holes in the halo, the dis- 
tribution should approximate the simple double Gaussian 
proposed by K. Schwarzschild (s ee BT08 jj4.4.3). The ver- 
tical velocity distribution (.Nordstrom et al\\2004 ). on the 
other hand, does have a re laxed appearance, as noted by 
ISeabroke fc Gilmord (I2007D . 

Various dynamical agents have been proposed to ac- 
count fo r kine matic features in the solar neighborhood. 
iKalnaisI ()1991h argued that the OLR of the bar in the 
MW might be close to the solar circle. Features in the 
subsequently-released HIPPARCOS data have also been at- 
tributed to the OLR of the bar ( Raboud et al. 1998; Dchncn 
l2000l : lFia3[200lh . while ISellwoodI (|20ia) attributed another 
to a recent ILR. 

De Simone et al. (|2004[ ) were able to produce distribu- 
tions of stars with a similar degree of substructure in sim- 
ulations of test particles moving in the adopted potential 
perturbations represen ting a succession of s hort-h ved spi- 
ral transients; see also lMinchev fc QuillenI ()2006f ). Other 
m odels tha t inclu d ed both bars and sp irals were presented 
bvlQuillenI (2003'), 'Chakrabarty' (2007), and'A ntoia et al\ 
(|2009D . while Hclmi et al. (2006) suggested that some of 
the substructure may also be caused by satellite infall (see 
also fTT|). 

Since spirals (and bars) are inefficient at exciting vertical 
motions, another scattering process needs to be invoked to 
redirect in-plane random motion into vertical motion. The 
influence of clumps in the disk, which are known to exist, 
seems consistent with the relaxed appearance of the verti- 
cal velocity distribution. Furthermore, the observed axes 
of the velocity ellipsoid, except possibly for the youngest 
stars (Fig. [16]), are consistent with the prediction of cloud 
scattering. Thus cloud scattering seems to be just suffi- 
cient to redirect velocities, which they are good at, but 
not to contribute significantly to heating. 



11. FRAGILITY OF DISKS 

Most of the stars in spiral galaxies reside in remarkably 
thin disks. In a photometric study of edge-on galaxies, 
lYoachim fc DalcantonI (|2006D found that the fraction of 
the baryons in the thin disk component is in the range 
70% to 90%, with higher fractions being characteristic of 
more massive galaxies. The so-called superthin disks are 
even more extreme; the disk of the low-luminosity galaxy 
UGC 7321, for example, has a radial sca le-length some 1 4 
times that of its characteristic thickness (jMatthewsll2000( ) . 

Hierarchical galaxy formation scenarios predict that galaxy 
formation is far from monolithic, with occasional major 
mergers with other halos, and frequent minor mergers. 
Stellar disks that have formed in t he centers of the ha- 
los ar e torn apart in major mergers ([Barnes fc Hernguistl 
119921 ). but the consequences of mino r mergers are harde r 
to determine from sim ulations (e.g. [Walker et al]ll996f ). 
[Toth fc Ostrikeil ([19921 ) argued that the existence of thin 
disks in galaxies can be used to constrain the rate of minor 
mergers. Their numerical estimates have been criticized on 
va rious grounds (e.g. Huang & Carlbcrg 1997; ScUwood et 
g/. [19981 : [Velazquez fc Whiteill999l ). but it is clear that a 
tight constraint remains. 

Wyse (this volume) stresses that the thick disk of the 
MW contains only very old stars, and that the ages of thin 
disk stars stretch back 10 Gyr. This fact would seem to 



imply that the last galactic merger to have stirred the MW 
disk occurred some 10 Gyr ago, and that no comparable 
disturbance could have occurred since. The MW may not 
be unique in this respect, as, Mould (2005,) found that thick 
disks in four nea r by ga laxies also appear to be old. 

Stewart et al. ([2008[ ) estimated the rate of mergers in 
the current ACDM cosmology. They concluded that 95% 
of parent halos of some 10^^ Mq will have merged 
with another halo of at least 20% of its mass in the last 
10 Gyr, and consequently the disk hosted by the parent 
must somehow survive in most cases. 

Possibly only a small fraction of infalling satellites pose 
a real threat to a disk. All satellites will be tidally stripped 
as they fall into the main halo, and some may even dissolve 
completely before they can affect the disk - th e Sagittar- 
ius dwarf galaxy appears to be a good example ([Law et al\ 
2005). A massive satellite loses orbital energy through dy- 
namical friction ( H9.8P causing it to settle deeper into the 
main halo. Its survival depends on its density (see BT08 
§8.3); it will be stripped of its loosely bound envelope until 
its mean density, p, is about one third the mean density of 
the halo interior to its orbital radius (BT08 eq. 8.92), and 
disrupted entirely as it reaches a radius in the main halo 
where even its central density falls below the mean interior 
ha lo density. This proces s has b een stu died in more deta il 
by'Boylan-Kolchin fc Ma (2007) and by Choi et all ([2009[ ). 
Thus the vulnerability of disks depends on the inner den- 
sities, in both dark matter and baryons, of the accreted 
sub-halos. 

If even a single, moderately massive core survives to the 
inner halo, it could cause an unacceptable increase in the 
disk thickness. Vertical heating of the disk occurs when 
a passing or penetrating satellite is able to increase the 
vertical motions of the disk stars. High speed passages 
therefore deposit little energy into disk motions, but if the 
satellite's orbit remains close to the disk mid-plane, then 
its vertical frequency will couple strongly to that of some 
of the disk stars and heating will be rapid. Indeed. Read et 
al. ([2008[ ) argued that the accretion event (s) that created 
the thick disk of the MW most probably resulted from the 
infall of a subhalo whose orbit plane was inclined at 10° 
- 20° to that of the disk. The energy deposited could 
take the form of exciting bending waves in the disk that 
can propag ate radially untfi they are damped at vertical 
resonances ([Sellwood aL[ [1998l). 

Kazai^ziHi^^ri 2009f ) simulated several minor merg- 
ers with sub-halos in detailed models with plausible pa- 
rameters. They found that the disk is substantially thick- 
ened and heated by the mergers, although they noted that 
their simulations lacked a gaseous component. Simula- 
tions that include a gas component cannot resolve the 
small-scale physical processes of gas fragmentation, star 
formation, feedback, etc., and therefore the behavior of 
the dissipative gas component necessarily includes some- 
what ad hoc prescript i ons f or these aspects of sub-grid 
physics. [Hopkins et al\ ([2009D stressed that gas in mergers 
settles quickly to begin to form a new disk; however, it 
is the fate of the stars that had been f ormed prior to the 
merge r that is the principal concern. [Kazantzidis et all 
([2009D argued it is possible that a dissipative component 
in the disk could absorb some of the orbital energy of the 
satellite which would reduce the heating of the stellar disk. 
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and such an effect appears to have occured in the simula- 
tions bv lMoster et al\ (|2010D . A clear conclusion has yet 
to emerge from this on- going research effort, but exclu- 
sively old thick disks together with the prevalence of thin 
disks poses a substantial, though perhaps surmountable, 
challenge to the ACDM model. 

It should be noted that ideas to reduce the density of the 
inner m ain halo thr ough frictional energy lo ss to halo {e.g. 
lEl-Zant et al. 20011: iMashchenko et aLll200Q: Romano-Dfaz 
et a/. l2008ar require the kinds of dense massive fragments 
that are themselves a danger to the disk. This will be of 
less importance in the early stages of galaxy assembly, but 
disk survival adds the requirement that any such process 
be completed quickly. 

12. CONCLUSIONS 

This review has focused on understanding the mecha- 
nisms that underlie the various instabilities and processes 
that affect the structure of galaxy disks. The discussion 
has referred extensively to simulations that are sufficiently 
simple that they clearly illustrate each particular aspect of 
the behavior. A complementary, and powerful, approach 
is to add additional physical processes to simulations, with 
the aim of improving realism to address otherwise inacces- 
sible questions. However, the increased complexity of the 
behavior that comes with increasing realism, makes a deep 
understanding of the results harder to achieve. Hopefully, 
the behavior described here will form a solid basis on which 
to build as the challenges presented by more realistic sys- 
tems are addressed. 

A contr ast with accretion di s k the ory seems appropri- 
ate here. IShakura fc SunvaevI ()1973f ) proposed a scaling 
relation for turbulent viscosity that led to rapid progress 
in modeling accretion disks some 25 ye ars before the likely 
origin of the viscosity was identified ()Balbus fc HawlevI 
|1998[ ). A theory for the structural evolution of galax- 
ies seems much farther away. While most of the impor- 
tant physical processes may have been identified, exactly 
how they drive evolution is still not fully understood, and 
even the evolutionary path remains vague. In the absence 
of this understanding, simplifying scaling laws cannot be 
identified with confidence and it seems best to work at 
improved understanding of the mechanisms at play. 

Galaxy dynamics has made immen se strides in the 45 
years since the chapter by lOortj (|1965( ) in the correspond- 
ing volume of the preceding series. Understanding of local 
wave mechanics and the mechanisms for the princi- 
pal global gravitational instabilities (§Sjll&ini) is well ad- 
vanced. Bending waves and global buckling modes (® 
are mostly understood, but not at quite the same level, 
while lop-sided modes {^E^ perhaps still require some more 
concerted effort. Understanding of bar structure O and 
the role of bars in galaxy evolution has developed beyond 
recognition, and disk heating f ijlOp seems more solidly un- 
derstood, with our confidence b eing strongly boo s ted b y 
the extremely valuable data from iNordstrdm et al\ (120041 1. 

But answers to a number of major questions of galaxy 
dynamics are still incomplete. Although the bar-forming 
instability is now understood (! j4.ip . it does not provide 
a clear picture of why only a little over half of all bright 
disk galaxies are barred f H9.2p . After bars, spirals arms 
are the second most prominent feature of disk galaxies 



and are probably responsible for the most dynamical evo- 
lution fifTO)). yet a deep understanding of their origin (^T]) 
remains elusive. The good progress made in recent years 
to understand the warps of galaxy disks (^H]) still has not 
supplied a satisfactory account of their incidence. 

Many of these outstanding issues may be bound up with 
how galaxies form, our understanding of which is currently 
making particularly rapid progress. 
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APPENDIX 

RELAXATION TIME IN SPHEROIDS AND DISKS 

A test particle moving at velocity v along a trajectory 
that passes a stationary field star of mass m with impact 
parameter b is deflected by the attraction of the field star. 
For a distant passage, it acquires a transverse velocity com- 
ponent \v±\ ~ 2Gm/(bv) to first order (BT08 eq. 1.30). 
Encounters at impact parameters small enough to produce 
defiections where this approximation fails badly are negli- 
gibly rare and relaxation is driven by the cumulative effect 
of many small deflections. 

If the density of field stars is n per unit volume, the 
test particle will encounter 6n = 2Trb6b nv stars per unit 
time with impact parameters between b and b + 6b. As- 
suming stars to have equal masses, each encounter at this 
impact parameter produces a randomly directed v± that 
will cause a mean square net deflection per unit time of 



, n f2GmY 

I ~ — — X 2'Kbob nv — 
\ bv J 



I-kG m?n 5b 
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The total rate of deflection from all encounters is the in- 
tegral over impact parameters, yielding 



db SnG'^m^n 



In A, (A2) 



where InA = ln(&max/&min) is the Coulomb logarithm. 
Typically one chooses the lower limit to be the impact 
parameter of a close encounter, femin ~ 2Gm/v^ , for which 
\v±\ is overestimated by the linear formula, while the up- 
per limit is, say, the half-mass (or effective) radius, R, of 
the stellar distribution beyond which the density decreases 
rapidly. The vagueness of these definitions is not of great 
significance to an estimate of the overall rate because we 
need only the logarithm of their ratio. The Coulomb loga- 
rithm implies equal contributions to the integrated deflec- 
tion rate from every decade in impact parameter simply 
because the diminishing gravitational influence of more 
distant stars is exactly balanced by their increasing num- 
bers. 

Note that the flrst order deflections that give rise to this 
steadily increasing random energy come at the expense of 
second order reductions in the for ward motion of the same 
particles that we have neglected (jHenonlll973f ). Thus the 
system does indeed conserve energy, as it must. 

We deflne the relELxation time to be the time needed 
for ~ u^, where v is the typical velocity of a star. Thus 



"^rela 



dnA 



(A3) 



To order of magnitude, a typical velocity ~ GNm/R, 
where iV is the number of stars each of mean mass m, 
yielding A ^ N. Defining the dynamical time to be Tdyn = 
R/v and setting N ^ R'^n, we have 



N 



F '^dyn 5 



(A4) 



61nA^ 

which shows that the coUisionless approximation is well 
satisfied in galaxies, which have 10® ^ N < 10^^ stars. 
Including the effect of a smooth dark matter component in 
this estimate would increase the typical velocity, v, thereby 
further lengthening the relaxation time. 
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This standard argument, however, assumed a pressure- 
sup ported quasi-spherical system in several places. Rv- 
bicki ()1972f ) pointed out that the flattened geometrv and 
organized streaming motion within disks affects the relax- 
ation time in two important ways. First, the assumption 
that the typical encounter velocity is comparable to the 
orbital speed v = {GNm/ RY/"^ is clearly wrong; stars 
move past each other at the typical random speeds in the 
disk, say jSv with /? 0.1, causing larger deflections and 
decreasing the relaxation time by a factor . 

Second, the distribution of scatterers is not uniform in 
3D, as was implicitly assumed in eq. (jAip . Assuming a ra- 
zor thin disk, changes the volume element from 2'kv bSb for 
3D to 2v Sb in 2D, which changes the integrand in eq. (IA2[) 
to b^^ and replaces the Coulomb logarithm by the factor 
(^min ~ ^max)- ^D therefore, relaxation is dominated 
by close encounters. Real galaxy disks are neither razor 
thin, nor spherical, so the spherical dependence applies at 
ranges up to the typical disk thickness, zq, beyond which 
the density of stars drops too quickly to make a signif- 
icant further contribution to the relaxation rate. Thus 
we should use A ~ zo/bmin for disks. More significantly, 
the local mass density is also higher, so that TV ~ R^z^n. 
These considerations shorten the relaxation time by the 
factor {zq/ R)h\{R/ zq). An additional effect of flattened 
distribution of scatterers is to determine the shape of the 
equilibrium velocity ellipsoid, as discussed in iil0.3l 

A third consideration for disks is that the mass distri- 
bution is much less smooth than is the case in the bulk 
of pressure supported galaxies. A galaxy disk generally 
contains massive star clusters and giant molecular clouds 
whose influence on the relaxation rate turns out to be non- 
negligible (see ^lOp . 



